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Abstract 

This paper illustrates the relationship between 
Transmission Loading Relief (TLR) and Available 
Transfer Capability (ATC), and present possible 
approaches to reduce the need to invoke TLR procedures 
while improving the utilization of the transmission 
system. In open-access electricity markets the 
reservations of transmission services are based on the 
data resulting from the off-line computation of ATC. 
When ATC data is inaccurate or misused, security 
violations can emerge and cause operators and operation 
planners to invoke TLR procedures. In this paper ideas 
are demonstrated on 7-bus, 26-bus and IEEE118 bus 
system and results have been presented and analyzed 1 . 

Keywords: Hour-ahead ATC, off-line ATC, Transmission 
Loading Relief ( TLR ). 

1. Introduction 

The concept of competitive market rather than regulated 
ones has become prominent in the past few years. 
Economists and political analysts have promoted the idea 
that free markets can drive down costs and prices thus 
reducing inefficiencies in power production. This change 
in the climate of ideas has fostered regulators to initiate 
reforms to restructure the electricity industry to achieve 
better service, reliable operation, and competitive rates. 
Deregulation of the power industry was first initiated in 
United Kingdom, followed suit in Norway, Australia 
partly in India. 

The North American Electric Reliability Council (NERC) 
has developed an Operating Manual that reduces a 
multitude of issues associated with electric power system 
operations [l]-[3]. This manual includes Transmission 
Loading Relief (TLR) procedure for dealing with issues 
that focus primarily on the relationship between 
interchange transactions and security violations. Since 
interchange transactions are typically created through the 



This study has been implemented on Power world 
simulator at CVR College of Engineering. 



use of Available Transfer Capability (ATC) information, 
there should be a very close relationship between the 
need for TLR and the accuracy and use of ATC 
information. This paper examines two important aspects 
of ATC accuracy and use. The first considers the non- 
simultaneous nature of posted numbers. The second 
considers the changing of numbers between areas with 
connecting paths [11]. 

Quantifying the capabilities of a transmission system for 
interchange has been of interest for over 25 years. These 
earlier methods focused primarily on Simultaneous 
Interchange Capability (SIC) rather than wheeling 
capabilities. The application of SIC was more of an 
emergency capability in case local generation was lost 
and massive import was required. The terminology of 
Available Transfer Capability (ATC) was created to 
quantify open access to transmission facilities. As such, it 
is used to reserve transmission services needed for 
economic transactions between suppliers and consumers. 
The theoretical foundations and computational burdens of 
ATC are essentially identical to those of SIC. 

A large percentage of ATC calculations performed today 
utilize linear load flow techniques to approximately 
predict changes in line loading in response to transfers 
and line outages. In their simplest form, these linear 
methods rely exclusively on distributions factors [16]. 
The Power Transfer Distribution Factor (PTDF) is a 
“power” version of traditional linear circuit analysis with 
current division. The Line Outage Distribution Factor 
(LODF) also uses traditional linear circuit analysis to 
compute these “large-change” sensitivities. While some 
work has been done on including dynamic constraints on 
transfers, the core ATC computation considers static 
constraints only. Currently, the following procedure 
referred to as off-line ATC, is used to provide ATC 
information for commercial purposes: 

1. The security coordinator builds a regional model 
for the expected peak conditions. The model 
includes power flow base case, contingencies, 
monitored elements and directions. 
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2. The security coordinator computes firm and non- 
firm ATC based on the set of proposed directions 
on a monthly, weekly and daily basis, and posts 
the results on the Open Access Same-time 
Information System (OASIS). 

3. Marketers use ATC data to propose transactions 
and make reservations of transmission services. 

4. Control areas verify and debug the OASIS values 
for the day, including specific system changes and 
operational considerations. 

5. As the time for transaction implementation 
approaches, control areas perform security 
analysis to identify potential problems. 

The process has been implemented within a daily 
framework due to important technical reasons related to 
the fact that ATC is a massive and complicated 
computational task, and that making its results available 
demands a large data integration effort. The same reasons 
have left many ATC issues unsolved and set challenges 
regarding ATC procedures and computational methods. 

2. Background of ATC 

The Available Transfer Capability (ATC) is a measure of 
the transfer capability remaining in the physical 
transmission network for further commercial activity 
over and above previously committed uses. 
Mathematically, ATC is defined as the Total Transfer 
Capability (TTC) less the Transmission Reliability 
Margin (TRM) and the Capacity Benefit Margin (CBM) 
[ 1 ]. 

. ATC = TTC _ XRM _ CBM 

TTC is defined as the amount of electric power that can 
be transferred over the interconnected transmission 
network in a reliable manner while meeting all of a 
specific set of defined pre and post-contingency system 
conditions. 

TRM is defined as that amount of transmission transfer 
capability necessary to ensure that the interconnected 
transmission network is secure under a reasonable range 
of uncertainties in system conditions. 

CBM is defined as that amount of transmission transfer 
capability reserved by load serving entities to ensure 
access to generation from interconnected systems to meet 
generation reliability requirements. 

References [7]- [9] & [12] provide additional analysis of 
errors associated with ATC. 

3. Transmission Loading Relief 

Transmission Loading Relief (TLR) is a sequence of 
actions taken during operations planning or during real- 
time operation to avoid or remedy security violations 
associated with the transmission system. 

The process starts when a security coordinator identifies 
a transmission facility within its security area that is 
about to, or has exceeded the operating security limit. At 
this point the security coordinator may invoke TLR. This 
NERC TLR procedure involves following “levels” [1]: 
Level 1: Notify reliability coordinators of potential 
operating security limit Violations. 

Level 2: Hold transfers at present level to prevent 
operating security limit Violations. 




Level 3a: Reallocation of Transmission Service by 

curtailing interchange transactions using non-firm point- 

point transmission service to allow interchange 

transactions using higher priority transmission service. 

Level 3b: Curtail interchange transactions using Non- 

Firm Transmission Service Arrangements to mitigate an 

operating security limit Violation 

TLR Level 4: Reconfigure Transmission 

TLR Level 5a: Reallocation of Transmission Service by 

curtailing interchange transactions using Firm Point-to- 

Point Transmission Service. 

TLR Level 5b: Curtail interchange transactions using 
Firm Point-to-Point Transmission Service to mitigate an 
operating security limits Violation. 

TLR Level 6: Emergency Procedures. 

TLR Level 0: TLR concluded. 

4. Types of Transactions 

In open-access electricity markets transactions can be 
done in three different ways. They are: 

A. Simultaneous transactions 

B. Sequential transactions: off-line ATC 

C. Sequential transactions-updated ATC 

A. Simultaneous transactions 

In a system, suppose that the transmission reservations 
occur simultaneously [4], [5]. If these transactions are 
implemented, there is a chance that one or more lines 
may be overloaded. When this schedule is analyzed prior 
to implementation, TLR will have to be invoked to 
relieve the overloading of these lines. Each participating 
control area for the transaction will implement a local 
TLR procedure, which would imply the rearrangement of 
its own resources, or ask the security coordinator for a 
TLR procedure, which would require the rearrangement 
of the resources of all the interconnection entities 
involved in the transaction. 

Once TLR is initiated, the process determines for each 
type of transaction the impact of those that affect the 
loading of the facility in 5% or more based on PTDFs [6]. 
Then transactions are curtailed in a progressive, 
prioritized manner. If no more transactions can be 
curtailed, the process goes to the next TLR level. 

B. Sequential transactions: off-line ATC 

As an alternative approach, consider a sequence of 
transactions starting from base case. In this case, ATC 
values remain as those computed by the off-line ATC 
process as the transactions take place. 

1) At time tl marketer Ml proposes a transaction 
in one direction. The transfer should create no 
problem if there is ATC available after the 
transaction is implemented. 

2) At time t2 marketer M2 reads ATC for another 
direction and proposes a transaction T2. There 
may be a chance of overload if the transaction 
T2 is made considering the ATC values prior to 
transaction Tl that needs to be relieved. So, 
TLR procedure should be invoked to relieve this 
overload. 

Note in this case that for proposed transactions the values 
of unreserved ATC do not change and therefore the 
capability for each direction is assumed to be the one that 
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was computed for the base case. As this is a non- 
simultaneous process, new transactions can overload the 
system. TLR will be needed even though the transactions 
were based on posted ATC values. 

C. Sequential transactions -updated ATC 

Because of the disadvantage of invoking TLR procedures 
in the Sequential Off-Line ATC method, we go for this 
case. Start again with the system at the initial operating 
point. 

1) At time tl marketer Ml proposes a transaction 
in one direction. The transfer should create no 
problem if there is ATC available after the 
transaction is implemented. 

2) ATC values are recomputed and updated. 

3) At time t2 marketer M2 reads ATC for another 
direction and proposes a transaction T2. As the 
transaction is made based on the updated ATC 
values, there should be no overload. 

4) ATC values are recomputed and updated after 
each transaction. 

This case is an ideal case in which ATC is updated each 
time a transaction takes place, is proposed, or a 
reservation is made. Note that ATC values remain 
positive. For this ideal case, no TLR procedure is 
required [17]. 

If a continuous ATC computation is possible, and 
transactions are based on updated ATC data we can think 
of transactions as moving within a secure transactional 
space. 

5. Case Study 

A. 7-Bus System 

The computations were tested in 7-bus system with 3- 
areas. The system has been divided into three areas, 
namely; Area A, Area B and Area C. Area A includes 
buses 1, 2, 3, 4 and 5. Area B consists of bus 6 whereas 
Area C consists of bus 7. ATC values at initial operating 
are shown in Table 1 [13]. 



Table 1. 7-Bus case, ATC Data 



Transfer 

Areas 


Transfer 

Buses 


ATC 


Limiting 

Line 


(MW) 


A-B 


1-6 


54 


1-2 


2-6 


265 


2-6 


4-6 


281 


2-6 


A-C 


1-7 


56 


1-2 


2-7 


106 


2-5 


4-7 


152 


4-5 


B-C 


6-7 


111 


6-7 



Given this system data, suppose that the following 
transmission reservations occur, which are all possible 
after the initial posting: 

Direction 1-6:50MW 
Direction 2-6:200MW 
Direction 4-6:200MW 

If these transactions are implemented line 1-2 will have 
over load of 22% and line 2-6 will have an over load of 
95%. When this schedule is analyzed prior to 
implementation, TLR will have to be invoked to relieve 
the overloading of these lines. Sequential ATC data for 
this example is given in Table 2. If all the ATC numbers 
were computed after t3 the results would be as shown in 



Table 3. Sequential transactions-update ATC data for this 
case is given in Table 4. 



Table 2. Sequential ATC Data 



Transfer 


to 


tl 


t2 


t3 


Area 


Bus 


ATC 


Limit 


ATC 


Limit 


ATC 


Limit 


ATC 


Limit 


MW 


Line 


MW 


Line 


MW 


Line 


MW 


Line 




1-6 


54 


1-2 


4 


1-2 


4 


1-2 


4 


1-2 


A-B 


2-6 


265 


2-6 


265 


2-6 


65 


2-6 


65 


2-6 




4-6 


281 


2-6 


281 


2-6 


281 


2-6 


81 


2-6 




1-7 


56 


1-2 


56 


1-2 


56 


1-2 


56 


1-2 


A-C 


2-7 


106 


2-5 


106 


2-5 


106 


2-5 


106 


2-5 




4-7 


152 


4-5 


152 


4-5 


152 


4-5 


152 


4-5 


B-C 


6-7 


111 


6-7 


111 


6-7 


111 


6-7 


111 


6-7 



Table 3. ATC for Final Conditions 



Transfer 


t3 


Area 


Bus 


ATC 


Limit 


MW 


Line 




1-6 


-168 


1-2 


A-B 


2-6 


-166 


2-6 




4-6 


-186 


2-6 




1-7 


-448 


1-2 


A-C 


2-7 


-438 


2-5 




4-7 


-513 


4-5 


B-C 


6-7 


87 


6-7 



If a continuous ATC computation is possible, and 
transactions are based on updated ATC data we can think 
of transactions as moving within a secure transactional 
space as shown in Figure. 1. 




Figure 1. Secure Transactional Space for transactions of 7-bus System 

Secure transactional space A secure transactional space 
(as mentioned earlier) is now created based on the above 
Sequential Updated ATC data. This space gives the limits 
for future transactions and all transactions are to be made 
within its limits to reduce the chance of overloading. 

It can be seen that each time a transaction is proposed, 
the secure transactional space changes (it can be reduced 
or augmented). 




3 




ACSE Journal, ISSN 1 687-481 1 , Volume 1 2, Issue 1 , Delaware, USA, June 201 2 



Table 4. Sequential Transactions-Update ATC 



Transfer 


to 


tl 


t2 


t3 


t4 


Area 


Bus 


ATC 


Limit 


ATC 


Limit 


ATC 


Limit 


ATC 


Limit 


ATC 


Limit 


MW 


Line 


MW 


Line 


MW 


Line 


MW 


Line 


MW 


Line 




1-6 


54 


1-2 


4 


1-2 


9 


1-2 


4 


1-2 


54 


1-2 


A-B 


2-6 


265 


2-6 


220 


2-6 


25 


2-6 


13 


2-6 


59 


2-6 




4-6 


281 


2-6 


28 


1-2 


26 


2-6 


14 


2-6 


62 


2-6 




1-7 


56 


1-2 


6 


1-2 


9 


1-2 


5 


1-2 


56 


1-2 


A-C 


2-7 


106 


2-5 


98 


2-5 


54 


2-5 


34 


2-6 


63 


2-5 




4-7 


152 


4-5 


34 


1-2 


57 


1-2 


30 


1-2 


90 


2-5 


B-C 


6-7 


111 


6-7 


123 


6-7 


68 


2-5 


69 


2-5 


79 


2-5 



For the transactional space shown in the above Fig. 1, the 
coordinates are the transactions in directions 1-6, 2-6 and 
4-6. The center of coordinates (0, 0, 0) represents the 
initial operating point where no transactions take place. 
At this point the secure transactional space is given by 
cube tO with dimensions {54, 265, 281}. Then T1 results 
in a secure transactional space given by cube tl with 
dimensions {4, 220, 28}. T2 moves the operating point to 
(50, 200, 0) and ATC results in a secure transactional 
space given by cube t2 with dimensions {9, 25, 26}. T3 
moves the operating point to (50, 200, 10) and so on. As 
it can be seen, each time a transaction is proposed, the 
secure transactional space changes (it can be reduced or 
augmented) after ATC values and the distances to each 
boundary are given by the ATC value for that direction. 

B. 26-bus System 

A 26-bus power system, taken from an example in a book 
by Sadat [10] is considered as a single area. Suppose that 
the following transmission reservations occur, which are 
all possible after the initial posting: 

Direction 1-18:200MW 
Direction 2-12:150MW 
Direction 4-9:150MW 

If these transactions are implemented, line 1-18 will have 
over load of 7% and line 7-9 and 4-12 will have an over 
load of 8%. When this schedule is analyzed prior to 
implementation, TLR will have to be invoked to relieve 
the overloading of these lines. Sequential transactions- 
update ATC data for this case is given in Table 5. Secure 
transactional space for this case is shown in Figure. 2. 

C. IEEE 118- Bus System 

For the case study, a large sized system, IEEE 118-bus 
system is composed of 54 generators with a total installed 
capacity of 8190MW and the system demand is 3668MW. 
The system is divided into four areas as follows. For this 
analysis each area is considered as one generating 
company. The details of areas are given in Tables 6. 

Table 5. Details of Areas of IEEE 118-bus S ystem 



Area 


Buses 


A 


1 - 38, 113 - 115 & 117 


B 


39-68 


C 


69- 102, 116 & 118 


D 


103-112 



Suppose that the following transmission reservations 
occur, which are all possible after the initial posting: 
Direction 25-59:700MW 
Direction 100-45:700MW 
Direction 80-ll:700MW 



If these transactions are implemented, line 8-5 will have 
over load of 24% and line 68-81 and 80-81 will have an 
over load of 31%. When this schedule is analyzed prior 
to implementation, TLR will have to be invoked to 
relieve the overloading of these lines. Sequential 
transactions-update ATC data for this case is given in 
Table 6. Secure transactional space for this case is 
shown in Figure. 3. 




Figure 2. Secure Transactional Space for transactions of 26-bus 
System. 




Figure 3. Secure Transactional Space for transactions of 118-bus 
System 
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Table 6. Sequential Transactions-Update ATC for 26-bus System 



Transfer 


to 


tl 


t2 


t3 


t4 


Bus 


ATC 


Limit 

Line 


ATC 


Limit 

Line 


ATC 


Limit 

Line 


ATC 


Limit 

Line 


ATC 


Limit 

Line 


MW 


MW 


MW 


MW 


MW 


1-18 


231 


1-18 


31 


1-18 


5 


1-18 


1 


1-18 


205 


1-18 


1-12 


240 


2-8 


74 


1-18 


15 


1-18 


3 


1-18 


165 


2-8 


2-12 


220 


2-8 


99 


1-18 


19 


1-18 


4 


1-18 


151 


2-8 


3-9 


113 


3-13 


93 


3-13 


40 


1-18 


7 


1-18 


86 


3-13 


4-9 


175 


12-10 


162 


10-12 


146 


1-18 


46 


1-18 


95 


7-9 


4-18 


190 


16-17 


42 


1-18 


9 


1-18 


2 


1-18 


112 


4-12 


26-18 


130 


5-6 


47 


1-18 


10 


1-18 


2 


1-18 


83 


11-26 



Table 7. Sequential Transactions-Update ATC for IEEE 118-bus System 



Transfer 


to 


tl 


t2 


t3 


Area 


Bus 


ATC 


Limit 

Line 


ATC 


Limit 

Line 


ATC 


Limit 

Line 


ATC 


Limit 

Line 


MW 


MW 


MW 


MW 


A-B 


25-59 


765 


64-65 


068 


64-65 


065 


64-65 


060 


64-65 


26-62 


841 


64-65 


352 


64-65 


335 


64-65 


246 


64-65 


A-C 


12-77 


794 


65-68 


862 


65-68 


965 


65-68 


1048 


65-68 


26-90 


836 


65-68 


909 


65-68 


1000 


89-90 


1007 


89-90 


B-A 


54-11 


861 


38-65 


835 


38-65 


770 


38-65 


665 


38-65 


65-27 


842 


38-65 


816 


38-65 


754 


38-65 


650 


38-65 


B-C 


54-77 


697 


65-68 


758 


65-68 


848 


54-56 


864 


68-69 


65-82 


580 


65-68 


630 


65-68 


705 


65-68 


673 


68-69 


C-A 


80-11 


801 


65-68 


739 


65-68 


90 


81-80 


89 


81-80 


92-32 


910 


65-68 


835 


65-68 


358 


81-80 


100 


81-80 


C-B 


87-54 


300 


86-87 


300 


86-87 


300 


86-87 


93 


81-80 


100-45 


808 


86-87 


700 


65-68 


327 


81-80 


77 


81-80 



6. Chaining ATC Numbers 

When a transaction involves several control areas and 
associated ATC postings, the need to reserve 
transmission services along a contract path can result in a 
misuse of ATC numbers. For example, if transmission 
services are required between areas A, B and C for an A 
to C transaction, reservations might be made based on the 
minimum capability posted for A to B and B to C. This is 
a logical assumption, but can result in a TLR procedure 
because this chaining is not normally valid. One of the 
primary reasons why it is not valid is because ATC 
numbers are based on maximum allowable transfers 
subject to numerous contingencies. It is very likely that 
the limiting transfer will occur because of a proposed 
contingency. Since the A to B transfer might be limited 
by one contingency and the B to C transfer might be 
limited by a different contingency, it is not correct to use 
the minimum of these two ATC numbers to determine 
the A to C maximum transfer. It may well be that the A 
to C transfer is actually limited by a completely different 
contingency than the A to B and B to C. Furthermore, the 
ATC numbers may be the result of different limiting 
elements. There is one case where the ATC number for A 
to C can be obtained from the ATC numbers for A to B 
and B to C, but the result is not simply the minimum of 
the two. Consider the following new line ij flows in 
response to a transaction of P from A to B, B to C and A 



to C: 






p new 
r ij, AB 


_ p old 

r ij 


. rp new 
1 ij, AB 


p new 
r ij, BC 


p old 

ij 


. rp new 
+ A ij, BC 


p new 
r ij, AC 


_ p old 

ij 


. rp new 
1 ij, AC 



For a maximum allowable line ij flow of Pij max , the 
maximum P which can be used for each transaction is: 



P 

max AB 



^max BC 



p 

A max AC 



(pma. _ p old } 
T ij,AB 

(P“ - Pif d ) 

T ij,BC 

(Pi“ - Pf ) 

Tij, AC 



Where all of the P’s and T’s are assumed to be positive 
numbers. In this case, P maX Ac is the correct ATC for a 
transaction from A to C. If only the ATC postings for A 
to B and B to C are given, the ATC for A to C can be 
computed from the other two as (by inspection from 
above, since 



ij, AC 



- T -- 



+ T 



ij, BC 



P 



max AB 



(i^max AB ^max BC ) 
(i^max AB ^max BC ) 



Note that this is not normally the same as the minimum 
of 

PmaxAB and PmaxBC- Also note that this computation 
may need to be revised when one or more of the T’s are 
negative numbers, or a line flow maximum is reached 
through a negative flow. Also note that this is for the 
special case where all ATC numbers are determined by 
the same limiting line and the same contingency. 

Further investigation into the use of posted ATC for 
chaining needs to be done. It may be that numbers other 
that ATC should be used. For example, perhaps PTDF’s 
should be posted rather that the ATC. Then ATC can be 
computed when needed for any transaction [15]. 

As a compromise between offline ATC computation and 
real-time ATC, the implementation of hour-ahead ATC 
would be consistent with hourly market decisions. This 
would solve a large number of the problems. The 
following section gives a framework for this. 
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7. Hour- Ahead ATC 

Hour-ahead ATC performs the same computation as off- 
line ATC. It operates though over the planned model for 
the next few hours. It would be used intensively within 
the Energy Management System (EMS) study mode, 
imposing hard performance and data integration 
challenges. 

Current ATC solvers are able to provide accurate results 
for systems with thousands of buses, and hundreds of 
directions and contingencies in few minutes. The difficult 
part of hour-ahead ATC is data management and 
modeling. 

To accurately compute ATC, the following information 
must be provided to the solver: 

System and Transaction data: network topology, 

equipment parameters, control settings, load profile, 
generation profile, external model description, 
interchange schedule, contingency description, 
subsystem description, direction description. 

From existing technologies, ATC for the next hour can 
be computed as follows. 

1) Obtain a base system case from the real-time 
state estimator solution. This would give a base 
topology description, generation schedule and 
base power flows. 

2) Include in the system model modifications to 
loads using system and bus load forecast 
functions. Update generation settings using 
generation schedules and update base flows 
using information for schedules and 
transmission services. 

3) Solve the power flow for the next hour case and 
compute PTDFs and LODFs matrices. 

4) Compute ATC for confirmed reservations using 
a transaction analyzer based on linear methods 
which could also be used for other computations 
such as contingency analysis, sensitivity 
analysis, interchange distribution calculations, 
transaction arrangement and transmission 
loading relief [14]. 

Hour-ahead ATC would run cyclically (every hour) over 
future hour models. It could also run on operator request 
and when new reservations have been confirmed [17]. 
The function would be the core of the information system, 
serving both the economic and reliability goals of the 
hourly-market. A flow chart of the functions is shown in 
Figure 4. 

8. Conclusion 

TFR procedures are implemented to relieve overloading 
in the transmission system. TFR is required because 
proposed transactions are based on ATC values that do 
not necessarily reflect the actual transfer capabilities 
either due to the non- simultaneous nature of ATC, or 
improper chaining. In order achieve an efficient and 
stable performance of the hourly market; the participants 
need to have accurate ATC information. A way to 
increase the quality and meaningfulness of ATC data is 
to compute it in such a time framework that the model 
included in the computation is the closest to the real 
situation. In the ideal case, marketers would have 
updated information each time a reservation is confirmed. 



With existing technologies it may be possible to 
incorporate hour-ahead ATC functions that minimize the 
need for TFR as well as enhance the utilization of 
transmission resources. The reduction of the need for 
TFR may also be achieved by minimizing misuse of 
ATC numbers, or through alternative concepts that are 
more accurate and flexible for use in reserving contract 
paths. 




I I Function 



Figure 4. Hour-Ahead ATC 
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Nomenclature 

Ty,AB Transaction of Real power for line ij from area A 
to area B 

Tij,B C Transaction of Real power for line ij from area B 
to area C 

Ty,Ac Transaction of Real power for line ij from area A 
to area C 

Pi j; old Real power in line ij before the transaction Ty 

Pij,AB new Real power in line ij from area A to B after the 
transaction T^ab 

Pij,Bc new Real power in line ij from area B to C after the 
transaction T^bc 

Pij,Ac new Real power in line ij from area A to C after the 
transaction T^ac 

p.. max Maximum allowable power flow in line ij 

PmaxAB Maximum power that can be transferred from area A 
to area B 

P maxBC Maximum power that can be transferred from area B 

to area C 

P max ac Maximum power that can be transferred from area A 
to area C 




6 




ACSE Journal, ISSN 1 687-481 1 , Volume 1 2, Issue 1 , Delaware, USA, June 201 2 



References 

[1] NERC Operating Manual (latest edition), 
available at http:www.nerc.com 

[2] Transmission Transfer Capability Task Force, 
Transmission Transfer Capability, North 
American Electric Reliability Council, 
Princeton, New Jersey, May 1995. 

[3] Transmission Transfer Capability Task Force, 
Available Transfer Capability Definitions and 
Determination, North American Electric 
Reliability Council, Princeton, New Jersey, June 
1996. 

[4] G. L. Landgren, H. L. Terhune, R. K. Angel, 
Transmission Interchange Capability - Analysis 
by Computer, IEEE Transactions on Power 
Apparatus and Systems, Vol. PAS-91(No. 6): 
2405-2414, Nov/Dec 1972. 

[5] G. L. Landgren, S. W. Anderson, Simultaneous 
Power Interchange Capability Analysis, IEEE 
Transactions on Power Apparatus and Systems, 
Vol. PAS-92(No. 6): 1973-1986, Nov/Dec 1973. 

[6] P. W. Sauer, On the Formulation of Power 
Distribution Factors for Linear Load Flow 
Methods, IEEE Transactions on Power 
Apparatus and Systems, Vol. PAS-100(No. 2): 
764-770, February 1981. 

[7] P. W. Sauer, Technical Challenges of 
Computing Available Transfer Capability 
(ATC) in Electric Power Systems, Proceedings 
of the Thirtieth Annual (1997) Hawaii 
International Conference on System Sciences, 
Vol. V, pp. 589-593, Maui, Hawaii, January 7- 
10, 1997. 

[8] P. W. Sauer, Alternatives for calculating 
Transmission Reliability Margin (TRM) in 
Available Transfer Capability (ATC), 
Proceedings of the Thirty-First Annual (1998) 
Hawaii International Conference on System 
Sciences, Vol. Ill, Kona, Hawaii, January 6-9, 

1998, p. 89. 

[9] S. Grijalva and P. W. Sauer, Reactive Power 
Considerations in Linear ATC Computation, 
Proceedings of the Thirty-Second Annual (1999) 
Hawaii International Conference on System 
Sciences, pp. 1- 11, Maui, Hawaii, January 5-8, 

1999. 

[10] H. Saadat, Power System Analysis, McGraw- 
Hill, Boston, 1999. 

[11] Masoud Rashidinejad, Hamid Khorasani, Amir 
Rashidinejad, Transmission Expansion Planning 
in Restructured Electricity Industry Using a 
Hybrid Heuristic Technique, International 
Review of Modelling and Simulations , Vol. 3 
(Issue 3); 283-289, June 2010. 

[12] S. Grijalva, P. W. Sauer and James D Weber, 
"Enhancement of Linear ATC Calculations by 
the Incorporation of Reactive Power Flows", 
IEEE Transactions on Power Systems, Vol. 18, 
No. 2, MAY 2003, pp. 619-624. 

[13] Power World Simulator. Version 8.0 Education 
Edition. Power World Corporation, Urbana, IL, 
USA. 



[14] Niimura, T., Niu, Y., “Transmission Congestion 
relief by economic load Management”, IEEE 
Power Engineering Society Summer Meeting, 
Vol. 3, pp. 1645-1649, July 2002. 

[15] Azhar B. Khairuddin, S. Shahawaz Ahmed, M. 
Wazir Mustafa, Abdullah A. Mohd. Zin and 
Hussein Ahmad, “A Novel Method for ATC 
Computations in a Large-Scale Power System”, 
IEEE Transaction on Power Systems, Vol. 19, 
No. 2, pp. 1150-1158, May 2004. 

[16] Paramasivam Venkatesh, Ramachandran 

Gnanadass and Narayana Prasad Padhy, 
“Available Transfer Capability Determination 
using Power Transfer Distribution Factors”, 
International Journal of Emerging Electric 
Power Systems, Vol. 1, Issue 2, Article 1009, 
2004. 

[17] S B Yokoyoma, R B Niimura and T A Niioka, 
“Transmission Loading Relief Solutions for 
Congestion Management”, Electric Power 
Systems Research, Vol. 67, Issue 2, pp. 73-78, 
February 2003. 



Ch Sunil Kumar born in khammam 
on February 25, 1977 received B. 
Tech degree in electrical and 
electronics engineering in the year 
2000 and M. Tech degree in 
Electrical Power Engineering the 
year 2006 from JNTU College of 
Engineering, Kukatpally, 

Hyderabad, AP, India. He is a PhD student at JNTUH in 
the Electrical & Electronics Engineering. His area of 
interest includes Electrical Power systems, Power 
electronics and Electrical Machines. Mr. Chava is life 
member of Indian Society for Technical Education 
(ISTE). 

Dr. J. Amarnath obtained the B.E 
degree in electrical engineering 
from Osmania University, 
Hyderabad, A.P. India in 1982 and 
the M.E. degree in power systems 
from Andhra University, 
Visakhapatnam in 19 84. He 
worked in Tata Electric Company, Bombay during 1985- 
1986. In 1987 he was a Lecturer in Andhra University 
for a period of 2 years and then he joined in Nagarjuna 
University for a period of 4 years as Lecturer. In 1992 he 
joined JNTU College of Engineering, Kukatpally, 
Hyderabad. Presently he is professor and head of the 
department of Electrical and Electronics engineering 
department, JNTU, Hyderabad, A.P. He presented more 
than 250 research papers in national and international 
conferences. His research interests includes high voltage 
engineering, gas insulated substations, industrial drives, 
power electronics, power systems, microprocessors and 
microcontroller applications to power systems and 
industrial drives. Dr. Amarnath is life member of Indian 
Society for Technical Education (ISTE). 



Biographies 







ACSE Journal, ISSN 1 687-481 1 , Volume 1 2, Issue 1 , Delaware, USA, June 201 2 



Dr.P S Subrahmanyan received 
his bachelor of Engineering in 
Electrical & Electronics 
Engineering in the year 1960 & 
Master’s Degree in Electrical 
Power Systems from Jawaharlal 
Nehru Technological University 
in the year 1977. He received his 
PhD from IIT Madras in the year 1983. He published 
a number of papers in National and International 
Journals, Conferences, and several text books. 
Basically from Electrical Engineering discipline, he 
cross- migrated to the field of Computer Science and 
Engineering. His areas of interest include Fault 
analysis, six phase system & six phase induction 
motors. Dr. Pasupati sadasiva Subrahmanyam fellow 
of The Institution of Engineers (India), fellow of 
National Federation of Engineers, Senior Member of 
IEEE, Member of Computer Society of India, Member 
of Indian Society for Technical Education. 





ACSE Journal, ISSN 1 687-481 1 , Volume 1 2, Issue 1 , Delaware, USA, June 201 2 




1C!®®? 



www.icgst.com 



asm 



Voltage Instability and Voltage Collapse as Influenced by Cold Inrush 

Current 



Youssef A. Mobarak, and Mahmoud M. Hussein 

Electric Engineering Department, Faculty of Energy Engineering, South Valley University, Aswan, Egypt 

y_a_mubarak @ hotmail. com 



Abstract 

The occurrence of voltage instabilities or voltage 
collapses depend on the duration of the persistence of the 
fault and on the type of fault, some faults lead to voltage 
instabilities, others lead to voltage collapse. Evaluation of 
fault durations causing occurrence of voltage instabilities 
or collapse is the main goal of this paper. This paper 
searches for the effect of peak inrush current and its 
duration periods initiation of voltage instability, with 
lagging and leading load power factors at certain loads 
buses. Also, this paper searches about the fault duration 
required for voltage collapse. The influence of the cold 
inrush current on this fault duration is studied. Form 
several tedious trials, it is noticed that, the phenomena of 
voltage collapse can be occurred with some line opening 
in the studied system. In this paper, the power system 
dynamic simulation program is developed for dynamic 
analysis of voltage stability. 

Keywords: Power System Stability - Voltage Instability - 
Voltage Collapse - Cold Inrush Current. 

1. Introduction 

Voltage stability is currently one of the most important 
research areas in the field of electrical power systems. 
With the increased loading and exploitation of the power 
transmission system being also due to improved 
optimized operation the problem of voltage stability and 
voltage collapse attracts more and more attention [1-2]. 
Voltage collapse can take place in systems or subsystems 
and can appear quite abruptly which requires the 
improved continuous monitoring of the system state [1]. 
The problem of voltage collapse may be simply 
explained by an inability of the power system to supply 
the reactive power or by an excessive absorption of 
reactive power by the system itself. It is to be understood 
as a reactive problem and it is strongly affected by the 
load behavior [2-3]. The change in voltage is so rapid that 
voltage controls devices may riot take corrective actions 
rapidly enough to prevent cascading blackouts [3-5]. The 
continual increase in demand for electric power has 
forced utility companies to operate their systems closer to 
the limits of instability [6-7]. 



Voltage instability concerns voltage fluctuations around 
nominal values. These fluctuations are either periodic or 
non periodic [8-9]. The main factor causing voltage 
instability is the inability of the power system to meet the 
demand for service power. The heart of the problem is 
usually the voltage drop that occurs when active and 
reactive power flow through the series inductive 
reactances of the transmission network [10-13]. A 
criterion for voltage stability is that, at a given operating 
condition for every bus in the system, the bus voltage 
magnitude increases as the reactive power injection at the 
same bus is increased. A system is unstable if for at least 
one bus in the system, the bus voltage magnitude 
decreases as the reactive power injection at the same bus 
is increased. In other words, a system is voltage unstable 
if V-Q sensitivity is negative for at least one bus [14-15]. 
Progressive drop in bus voltage can also be associated 
with rotor angles going out of step. For example, the 
gradual loss of synchronism of machines as rotor angle 
between two groups of machines approach or exceed 1 80 
degree would result in very low voltages at intermediate 
points in the network. Voltage instability is essentially a 
local phenomenon; however, its consequences may have 
widespread impact [14]. Voltage collapse is a rapid 
progressive voltage fall and settling at certain value 
defined by system parameters, it is more complex than 
simple voltage instability and usually the result of a 
sequence of events accompany voltage instability leading 
to a low voltage profile in a significant part of the power 
system lasting long periods [16-17]. Features and 
modelling capabilities of conventional transient stability 
time domain simulation programs has been greatly 
enhanced over recent years to make them suitable for the 
assessment of voltage stability problems. 

In this paper, the power system dynamic simulation 
program is developed for dynamic analysis of voltage 
stability. The developed program has the modelling 
capability to account for dynamics, which are important 
in voltage stability analysis including those associated 
with dynamic loads and generators and their associated 
control. In this paper, studies the effect of the peak inrush 
currents as a disturbance on the duration period which 
causes voltage instability, with different lagging and 
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leading loads power factors. When inrush current is 
considered as a disturbance by increasing its peak value 
for many times of the normal operating system current, 
the duration period of this inrush current at which 
voltages will be unstable is determined for various values 
of initial peak inrush current from that network. The 
effect of the peak inrush currents on the fault duration 
causing to reach voltage collapse, when after certain 
steady state is reached, the line between two buses is 
opened, with different lagging and leading loads power 
factors. 

2. Power System Model Representation 

The power system considered consisting of n 
synchronous generators feeding through a transmission 
network, a number of loads. The system motion using 
one-axis generator model, under a disturbance, the 
following set of differential equations described for each 
generator [18]: 



Si =2 TTfoi CO . 


( 1 ) 




( 2 ) 




( 3 ) 



Where: S, f 0 , co, H, Z), and e . are the rotor angle, initial 

frequency, speed deviation, inertia constant, mechanical 
damping coefficient, and the q-axis voltage component, 
respectively. P m , P e , E fd , I d , and T d0 are the generator 

input mechanical power, output electrical power, 
excitation voltage, d-axis current component, and d-axis 
transient open-circuit time constant, respectively. Finally, 
X J and x\ are the d-axis and transient d-axis reactance’s. 

a a 

Note that, P e is computed as: 

P ei = E qi I qi + E di 1 di ( 4 ) 

Where, E d is the d-axis component of the voltage E , and 
/ is the q-axis component of the generator current, which 
is given as: 

/,. = (£■;- v)/(X) (5) 

Where, E is the voltage behind the reactance x d , and V 
is the generator terminal voltage. 

In this paper, we concerned with the study of the quasi 
static loads which it consists of heating and lighting 
equipments. Also, three quasi static models are used for 
induction motor load representation. The commonly used 
representations of static loads are either constant 
impedance to ground, constant current and constant real 
and reactive power. These load models are given in the 
following, where P L , Q L and V L are the load active power, 
reactive power, and bus voltage, respectively (their initial 
values are P Lo , Q Lo and V Lo ). 

Constant Impedance Model (C.Z.): 




Constant Current Model (C.I.): 




Constant Power Model (C.P.): 




Three models are used for induction motor load, and the 
active and reactive power has expressed by the 4 th order 
polynomials representation [14]: Induction motors with 
constant mechanical loads torque (T=Const.): 

P-0.275+8. 79V -23. 2 7V 2 +23.233 V 3 -8. 01 7V 4 
Q-- 16. 11+92.94V-179. 72V 2 +148. 56V 3 -44.67V* (9) 

Induction motors with mechanical load torque 

proportional to speed (T a co): 

P=-33+161. 67V -286. 033 V 2 +233. 033 V 3 -64.67V 4 
Q=48. 83-175. 67V+237.5V 2 -142.22V 3 +32. 56V* (10) 

Induction motors with mechanical load torque 

proportional to square of speed (T a a> 2 ): 

P-9. 233-81. 567V+309. 9V 2 -598. 5 V s +62 7. 967V 4 - 
342V 5 +75.96 7V 6 

Q=46. 739-53. 15V+ 172 1.67V 2 -325V 3 +3488. 704V 4 - 
1900V 5 +422.037 V 6 (11) 

The interrelation of the system elements is shown in 
Figure (1). The three systems solution is the dynamic 
simulation program, the network reduction program and 
the load flow program with loads represented by 
polynomial models. They should be used in sequence to 
solve the system equations. A flow chart of multi- 
machine power systems [49], written in MATLAB 

TOOLBOX is shown in Figure (2). 



Load 




Generator, 

Excitation, 

Turbine-Governor 



Figure 1 . Data requirements of power system simulation program 
and its support programs 



ACSE Journal, ISSN 1 687-481 1 , Volume 1 2, Issue 1 , Delaware, USA, June 201 2 



C Start j 



Read in: the parameters and operating 
conditions 



1 



Calculate initial 
conditions 

I 



Read in: Terminal hybrid matrix 

T 



Calculate load flow prior to disturbance 

— I — 





Compute the voltage 
for n loads to used it 
in the power equations 



Print results 



Solve the nonlinear equations for the m 
generators w.r.t references 




Figure 2. Digital simulation of multi-machine systems, flow chart 



The dynamic simulation program reads the nodal 
initial conditions, equivalent matrix, the generator, 
excitation, turbine-governor and finally the type of 
loads data. The output of the dynamic simulation 
program is the system open loop performance. The 
system performance is determined by solving the 
machine mathematical models together with the 
constraints imposed by the network. The non-linear 
machine models are solved numerically using the 
appropriate integration technique. The system 
performance is determined by solving the machine 
mathematical models together with the constraints 
imposed by the network. The non-linear machine 
models are solved numerically using the appropriate 
integration technique. 



For network having total N buses the calculated 
voltage at any bus is described by the following 
nonlinear equations: 

/ A 



1 



P K ~ JQk 



■Sw 



7=1 



( 12 ) 



With: K=l, 2, N, and K^S, where S means 

slack bus, £^is voltage at node K , £*is E K 
conjugate, y^Self admittance at node K , The 

mutual admittance between node K and node j, £ 
Voltage at node j, p k ,q k Scheduled real and reactive 
power entering the system at bus K. 

Pj i = P] co (CLo~^~CljS-\-Cl2S -\-Cl3S -\-U4S + ) 

Q k =Q k o(bo+b,S+b 2 S 2 +b 3 S 3 +b 4 S 4 + ) (13) 

Where: S=V/V 0 , and the active and reactive power 
components at the slack bus are then computed as: 



P S JQs ~ ^ ]PsK JQSK 



K= 1 
K*S 



(14) 






K= 1 
K*S 



Where: P s - jQ s Load at the slack bus, y sk Mutual 



admittance between S and K nodes, and Voltage 



at the slack bus. 

From all the dynamic models involved in this type 
of simulation, the electric alternators is the one 
requiring more calculations since the rest of the 
dynamic components are given by a block diagram, 
in the initial condition calculations are made by 
simply setting to zero any term containing a 
derivative term [14]. Using the phasor diagram 
shown in Figure (3) the initial condition for the 
synchronous generator can be calculated as follows: 
At the terminal side of every alternator, the known 
variables given by a load-flow study, are: Terminal 
voltage v , and Generated active power P, and 

* a 

reactive power Q. 




Figure 3. Phasor diagram of alternator initial condition. 
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Knowing these variables, it is easy to calculate the 
armature current of every unit and the corresponding 
power factor by: 



I a — 



f P + .iQ' 



(15) 



Where: / is the armature current phasor, and Resolve 
/ into components with V a as a reference: 
I r = I a cos (/), I x = —I a sin (j) ? and the phasor E qa in Figure 



(3), is given by: 



E qa ~V a +(r + jX q )l a 

= V a+ (^r + jh){ r + j X q) ( 16 ) 

= (V a -*qI X +rt r )+j(x q I r +rI x ) 



Then the angle (<5 - ft) is given by: 

arctan (I x r + I r x) 

o-p = y = — 

<V a -x q I x + rI r ) 



(17) 



where angle 5 is the position of the rotor. Then, calculate 
the terminal voltage d-q components: 

= -V a sin{S - p\V q = V a cos (S - fi) (18) 



Also, the armature current d-q components: 
h =- J a sin(£-/? + 4^ 

= I a cos(S-fi + </>) 



The field voltage of the machine from the stator side can 
be calculated by: E = V + rl - x,I, = E . The initial 

J q q CL CL rLf 

position of the rotor is calculated by: 8=y+P . 
Calculate the initial value of the state variable of the 
alternator model using the quantities just calculated. 
Repeat the whole process for all the electric alternators 
present. Once the bus admittance matrix representing the 
network has been built and its entries organized as 
generator nodes, non-linear load nodes and other nodes, it 
is necessary to eliminate those nodes with zero current 
injections. The procedure to obtain the reduced network 
is as follows: Assume the Y matrix is partitioned in 

the following way: 



V 




Y Y Y 

a gg ± gl ± GR 






h 


= 


Y Y Y 

*LG ± LL l LR 




v L 


0 




l 

1 




Vr_ 



Where: G means generator nodes, L means non-linear 
load nodes, R means remaining nodes. 



Eliminate the remaining nodes by successive elimination 
procedure: 



1 G 

h 



GG 



l LG 



l GL 



l LL 





1 

1 


1 


Vl_ 



Where: Y GG — Y GG 
Y 1 =Y - 

1 GL 1 GL 



-Y Y~ l Y 

1 GR 1 RR 1 RG 

Y Y~ l Y 

1 GR 1 RR 1 RL ’ 



Y 1 — y - Y Y~ l Y 

1 LG 1 LG 1 LR 1 RR 1 RG 

Y 1 =Y -Y Y~ l Y 

1 LL 1 LL 1 LR 1 RR 1 RL 



( 21 ) 



After adding the internal impedance of each generator to 
new admittance matrix, the equation (21) becomes re- 
arranged as follows: 




( 22 ) 



Where just y 1 change toy" , and eL means the voltage 

behind the transient reactance. The equation (22) can be 
re-arranged as follows, this is because the initial load 
current at load buses are known, and it is used to 
determine the load voltage at load buses respectively at 
the next step. 



X 




Xo 




l! 

G 




y L _ 




Hw 


Z LL_ 


h . 


(23) 



Where: Y* New admittance matrix, K,. r H T ^ Non- 

CjCj CjL ’ LCj 

dimensional matrices, Z LL Load impedance matrix. The 



equivalent matrix for the entire network is represented by 
the single line diagram of power system operating at the 
nominal loading condition is found as: 



Where: I, = — — ^ 



(24) 



PL = Pko( a 0+ a lSL+ a 2SL +CI3SL +CI4S1 4 + ), 

Gl=G^(^0+^7^L+^25 , L 2 + ^35'l 5 + ^A 4 + ) (25) 

The system differential equations are depicts in 
Appendix. 



3. Studied System 

The power system used for digital simulation consists of 
multi-machine nine-bus system developed by WSCC in 
United States. A single line impedance diagram of the 
system is shown in Figure (4), where the system is 
basically composed of three generating units and three 
loads, load A, load B, and load C are located at buses # 4, 
# 5, and # 6, respectively. Unit one is hydroelectric, while 
units two and three are steam driven generators. The 
system is considered operating normally for lsec. before 
fault is applied at # 5. The duration period which cause 
initiation of voltage instability in the considered system. 




12 




ACSE Journal, ISSN 1 687-481 1 , Volume 1 2, Issue 1 , Delaware, USA, June 201 2 




4. Cold Inrush Current 

Cold inrush current is described as the magnitude of 
instantaneous current drawn by the load when startup the 
load, it’s of short duration. Often milliseconds under a 
combination of certain conditions, when switching on 
solid state equipment such as computers, copiers, and 
electronic ballasts found in lighting system, as well as 
magnetic devices like motors, drives, core and coil 
ballasts, peak inrush current be several orders of 
magnitude greater than the circuit’s steady state current 
will be drawn. It can be as high as 20 times normal 
operating current. Usually this inrush current is reduced 
after 10 ms, but could take up to 30 to 40 cycles until the 
current is at a normal value. Switching power supplies 
produce high inrush current at turn on resulting from 
filter capacitor impedance. These large filter capacitors 
act like a short circuit, producing an immediate inrush 
current with a fast rise time [18-22]. 

High inrush current can affect on electrical components 
such as tripping of circuit breakers as well as arcing and 
failure of primary circuit components, such as switches. 
High inrush currents also necessitate over sizing of fuses 
or breakers which complicates other aspects of approves 
testing and may compromise protection on other vital 
components. Another side effect of high inrush is the 
injection of noise and distortion back into the mains. 
During startup, momentary connect bouncing in switches 
or relay may cause the contacts to become pitted due to 
arcing between the contact points. This surge in current 
can also cause serious damage, such as welding switching 
contacts together. High inrush current may shorten life of 
switches and motion sensors and other devices which are 
used break contact in the electrical system. It can cause 
heavy demand of reactive power and voltage problems at 



starting of network energizing. Starting of equipments 
and motors which require high starting currents than the 
normal values may suffer problems and failures. 

The different types of load and their special inrush 
current are the two important factors to the switching 
frequency of the relay, in addition, they also greatly 
affect the contact welding, particularly the inrush current 
should be taken into consideration, and a range that is 
safety enough should be granted. From the available 
literature [19], the different types of load and their 
corresponding inrush current are given in Table (1), and 
Figure (5) indicates some forms of inrush current for 
some types of loads. Inrush current can be evaluated at 
many times of rated load current. High inrush current is a 
times-dependent phenomenon caused by a coincidental 
set of circumstances that must occur simultaneously at 
the moment of switch on [22]. Measuring this high 
startup current is important to determine the proper 
corrective equipment to install, such as surge limiters, 
soft start devices or simply increasing wire size, special 
equipments are required. 

5. Cold Inrush Current Representation: 

The cold inrush currents power factors are taken as (0.5, 
0.6, 0.7, 0.8, and 0.9) pu. lagging or leading, and unity 
power factor. With each value of the load power factor, 
its recalculated. The inrush current has a value at the 
beginning and damp gradually to reach normal operating 
current. In our analysis, the inrush current decay takes the 
triangle form, as shown in Figure (6). Where: i p is the 
peak inrush current (i p — mi 0 ), i Q is the normal operating 
current, and t p is the period of decay of the inrush current 
in seconds. 
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Table (1) Loads inrush currents periods and their ratio’s to steady state values 



Type of load 


Inrush current / Rated current (i/i 0 ) 


Period of inrush current 


Resistive 


Steady state current 


0.07 to 0.1 second 


Solenoid 


10-20 times the steady state current 


0.07 to 0.1 second 


Motor 


5-10 times the steady state current 


0.2 to 0.5 second 


Incandescent lamp 


10-15 times the steady state current 


Approx. 1/3 second 


Mercury lamp 


3 times the steady state current 




Sodium vapor lamp 


1-3 times the steady state current 




Capacitive 


20-40 times the steady state current 


1/120 to 1/30 seconds 


Transformer 


5-10 times the steady state current 




Fluorescent lamp 


5-10 times the steady state current 


10 seconds or less 


Electromagnetic contact 


3 times the steady state current 


1/60 to 1/30 seconds 




(a) Incandescent lamp load 
(i/i o =10 -15 times) 



(b) Mercury lamp load 
(i/i 0 =3 times) 



(c) Fluorescent lamp load 
(i/i 0 =5 to 10 times) 




(d) Motor load (i/i 0 =5 to 10 times) 



(e) Solenoid load (i/i o =10 to 20 times) (f) Electronic contact load (i/i 0 =3 times) 



Figure 5. Some forms of inrush current for some types of loads 



Indicates the inrush current decay and its magnitude at 
any instantaneous through this decay till it reaches to 
steady state current by equation i-(t( i 0 -mi 0 )/t p ) +mi 0 , 
where: m is the peak inrush current times to its steady 
state value, t is the instantaneous time from ( 0 to t p ), and i 
is the instantaneous current. 




Figure 6. The form of inrush current in our analysis 

6. Considered Disturbances 

The considered disturbance is suddenly increasing of the 
peak inrush current many times of normal operating 
system current period of time at load bus # 5, with 
lagging and leading loads power factors at these buses 
successively, until the voltage instability starts to occur. 



The loads are represented by CZ loads at all load buses. 
The other disturbance, the time response for load buses 
voltages at load buses with varying the peak inrush 
current by many times of normal operating current such 
as (100%, 300%, and 500%) of normal operating current 
at load bus #5, at t=l sec., for a duration period equal 0.2 
sec., and opening the line connecting between two buses 
#5 and #7 at t=5 sec. until the voltage collapse happen, 
with lagging and leading loads power factors such as (0.8 
pu.). The system is let to reach steady-state operating 
condition opening at t=5 sec. The fault is taken as 
opening the line which connects between nodes #5 and 
#7 and lasts until state of voltage collapse initiation is 
attained. 

7. Results and Discussions 

7.1 Effect of suddenly increasing of the peak inrush 
current on the duration period for voltage instability: 

The duration period of inrush current which cause 
voltage instability when the inrush current suddenly 
increases, where it takes different initial peaks such as 
(16, 15, and 12) times of normal operating system 
current, at bus (5). Load power factors at load bus (5) is 
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taken as (0.5, 0.6, 0.7, 0.8, and 0.9) pu. lagging or 
leading. It can be noticed that, the duration period of 
suddenly increases of inrush current which cause voltage 
instability initiation increases by decreasing the value of 
the peak inrush current for each value of lagging, leading 
and unity load power factor at load bus (5), also it can be 
observed that, the duration period of suddenly increases 
of inrush current which cause voltage instability initiation 
increases by increasing the values of lagging power 
factors at load bus (5), and the duration period of 
suddenly increases of inrush current which cause 
initiation of voltage instability decreases by increasing 
the values of leading power factors at load bus (5) for the 
same value of peak inrush current, as shown in Figure 7 
(a, b, and c) for load power factors (0.5, 0.7, 0.9) pu. 
(lagging and leading), respectively. 

Also, Figure 8 (a, b) shows voltage stability limit, and 
indicate the zone of the unstable voltage and that of 
stable region with suddenly increasing of the peak inrush 
current by (12, 15, and 16) times of normal operating 
current system, with lagging and leading power factor 
respectively, at load bus (5) from 0.5 pu. to unity. Where 
the region of stability increases by decreasing the value 
of peak inrush current at any load power factor lagging or 
leading, and indicate that the duration period of suddenly 
increasing an inrush current which cause voltage 
instability initiation increases from lagging to leading 
loads power factor. 

7.2 Time response for loads voltages with different 
times of normal operating current system, due to 
suddenly increasing of (i/i 0 ) : 

The time response for load buses voltages at load buses 
with varying the suddenly increasing of the peak inrush 
current by many times of normal operating current 
system such as (16, and 12) times at load bus (5) as a 
disturbance, with lagging and leading load power factors 
such as (0.7 pu. Figure 9 indicates that, the time response 
for load buses voltages with suddenly increasing the 
value of peak inrush current of (16, and 12) times of 
normal operating current system, with 0.7 pu. lagging 
power factor at load bus (5). with the value of peak 
inrush current is 16 times of normal operating current, 
when the duration period is 200 msec., the system is 
stable and when the duration period is 201 msec., the 
system starts to be unstable, as shown in Figure 9(a, b), 
respectively. While, with the value of peak inrush current 
of 12 times the normal operating current, the duration 
period is 239 msec., and system voltage is stable, while 
when the duration period is 240 msec., the system voltage 
become unstable, as shown in Figure 9(c, d), 

respectively. 

Figure 10, we can observe that, the time response for 
load buses voltages with suddenly increasing the value of 
peak inrush current of (16, and 12) times of normal 
operating current, and with 0.7 pu. leading power factor 
at load bus (5), the duration period of suddenly increases 
of inrush current is 271 msec., and the system voltage is 
stable, while when the duration period of suddenly 
increases of inrush current is 272 msec., the system 
voltage starts to be unstable, as shown in Figure 10(a, b), 



respectively. While, with suddenly increasing the value 
of peak inrush current of 12 times of normal operating 
current, the duration period is 488 msec., and the system 
voltage is stable, while when the duration period is 489 
msec., the system voltage become unstable, as shown in 
Figure 10(c, d), respectively. Where: T c refers to duration 
period of inrush current. 

7.3 Voltage stability limits due to opening the line 
between with peak inrush current at load bus by 
(100%, 300%, and 500%) of normal operating 
current system: 

Figure 1 1 indicate the fault duration which cause 
initiation of voltage collapse when opening the line 
connecting between two buses (5) and (7) at t=5 sec., 
with the peak inrush current to normal load current is 
(100%, 200%, 300%, 400%, and 500%) for a period of 
time equal 0.2 sec. at t= 1 sec., at bus (5), load buses 
power factors are taken as (0.6, 0.7, 0.9) pu. lagging or 
leading power factors. All load buses are constant 
impedance loads. Also, the fault duration which cause 
voltage collapse decreases by increasing the peak inrush 
current for each value of lagging, leading and unity loads 
power factor at load bus (5), for the same peak inrush 
current. The fault duration which cause voltage collapse 
increases by increasing the lagging loads power factor 
towards unity power factor. But, it decreases slightly by 
increasing the leading loads power factor at load bus (5), 
as shown in Figure 1 1 (a, and b) for loads power factors 
(0.6, 0.7, 0.9) pu. (lagging and leading), respectively. 
Once more the decrease of the fault duration period is 
very slight with leading power factors. This is the same 
with higher lagging power factor loads. 

Figure 12 (a, b) shows voltage stability limit, and indicate 
the zone of the unstable voltage and that of stable region 
when opening the line between two buses (5) and (7) 
after the system reach to steady state operating condition, 
but in the beginning, the system suffer from increasing 
the peak inrush current of (100%, 300%, and 500%) of 
the normal operating current system can take a period as 
0.2 sec., with lagging and leading loads power factor, 
respectively at load bus (5) from 0.5 pu. to unity. Where 
the region of stability increases by decreasing the value 
of peak inrush current at any load power factor lagging or 
leading, and indicate that the fault duration period which 
cause voltage collapse initiation increases from lagging to 
leading loads power factor. 

7.4 Time response for loads voltages, with line 
opening, when (i/i G ) at load bus # 5, at 100%, 300%, 
and 500%: 

This section, shows that, the time response for load buses 
voltages at load buses with varying the peak inrush 
current by many times of normal operating current such 
as (100%, 300%, and 500%) of normal operating current 
at load bus (5), at t=l sec., for a duration period equal 0.2 
sec., and opening the line connecting between two buses 
(5) and (7) at t=5 sec. until the voltage collapse happen, 
with lagging and leading loads power factors such as (0.8 
pu.). Figure 13 indicates that, the time response for the 
load buses voltages with the peak inrush current is 
(100%, 300%, and 500%) of the normal operating 
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current, with 0.8 pu. lagging power factor at load bus (5), 
where, without inrush current (100% of the normal 
operating current), the fault duration which cause voltage 
collapse is 13.46 sec.. When peak inrush current is 
(300%), the fault duration which cause voltage collapse 
becomes 12.08 sec., and when the peak inrush current is 
(500%), the fault duration which cause voltage collapse 
goes to 10.92 sec., as shown in Figure 13(a, b and c), 
respectively. 



Figure 14 indicates the time responses load buses 
voltages with peak inrush currents (100%, 300%, and 
500%) of normal operating current, with 0.8 p.u. leading 
power factor at load bus (5), where, without inrush 
current (100% of normal operating current), the fault 
duration which cause voltage collapse is 22.02 sec., while 
when the peak inrush current is (300%) of normal 
operating current, the fault duration which cause voltage 
collapse becomes 21.16 sec.. For inrush current of 
(500%), the fault duration which cause voltage collapse 
turns to be 19.88 sec., as shown in Figure 14 (a, b and c), 
respectively. 






(c) power factor = 0.9 pu 

Figure 7. Effect of suddenly increasing of the peak inrush current many 
times at load bus (5) on the duration period for voltage instability, with 
lagging and leading power factor 




(b) with leading power factor 

Figure 8. Voltage stability limit, due to suddenly increasing of peak 
inrush current at load bus (5) by (12, 15, and 16) times of normal 
operating current system 




(a) i/i G = 16 times, and T c = 200 msec. 
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(d) i/i 0 = 12 times, andT c = 240 msec. 

Figure 9. Time response for loads voltages, with 0.7 pu. lagging pf., due 
to suddenly increasing of (i/i 0 ) at load bus (5) 



Figure 10. Time response for loads voltages, with 0.7 pu. leading pf., 
due to suddenly increasing of (i/i G ) at load bus (5) 




0 5 10 15 20 



Time (sec.) 




100% 200% 300% 400% 500% 

The inrush current (of normal operating current) 



(a) i/i G = 16 times, and T c = 271 msec. 



(a) with lagging power factor 
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(b) with leading power factor 

Figure 1 1 . Effect of peak inrush current wrt nominal current at load bus 
(5) on the fault duration for voltage collapse, when opening the line 
connecting between two buses (5) and (7), 




Lagging power factor at bus 5 (p.u.J 



(a) with lagging power factor 
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Figure 12. Voltage stability limit, due to opening the line between two 
nodes (5) and (7), with peak inrush current at load bus (5) by (100%, 
300%, and 500%) of normal operating current system a period 0.2 sec., 





(b) i/i 0 =300%, and T c =12,08 sec. 




(c) i/i 0 =500%, and T c =10.92 sec. 

Figure 13. Time response for loads voltages, with 0.8 pu. lagging power 
factor, with line between nodes (5) and (7) opening, when (i/i 0 ) at load 
bus (5) 




Time (sec.) 



(a) i/i 0 =100%, and T c =22.02 sec. 
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Time (sec.) 

(c) i/i o =500%, and T c =19.88 sec. 

Figure 14. Time response for loads voltages, with 0.8 pu. leading power 
factor, with line between nodes (5) and (7) opening, when (i/i 0 ) at load 
bus (5) 

8. Conclusions 

We can conclude that the voltage instability is influenced 
by cold inrush current application period: The duration 
period of increasing peak inrush current suddenly which 
cause voltage instability initiation increases by 
decreasing the value of peak inrush current, this is true 
with lagging and leading loads power factors. The 
duration period of increasing peak inrush current 
suddenly which cause initiation of voltage instability 
increases by increasing the loads power factor (in the 
lagging rang) for the same value of the peak inrush 
current. Voltage instability initiation decreases by 
increasing the leading loads power factor for the same 
value of the peak inrush current. 

The fault duration which cause initiation of voltage 
collapse is affected by the cold inrush current magnitudes, 
and the load power factors either in the leading or in the 
lagging ranges. For a sudden line opening, the fault 
duration period which cause voltage collapse initiation 
increases by decreasing the peak inrush current 
magnitudes for each value of lagging, leading and unity 
loads power factors. The fault duration which cause 
initiation of voltage collapse increases by increasing the 
lagging loads power factor, for the same value of the 
peak inrush current. While, the fault duration period 
which cause voltage collapse decreases by increasing the 
value of leading loads power factor for the same value of 
the peak inrush current magnitudes. 
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Abstract 

This paper proposes and validates models to 
accurately represent Static Synchronous 
Compensators STATCOM regulates system voltage 
in transient (load injection and rejection) by 
absorbing or generating reactive power of sample 
three-bus test system. The proposed STATCOM 
models are justified based on the basic operational 
characteristics of this Flexible AC Transmission 
System FACTS controller for both phase and PWM 
control strategies. The complete digital simulation of 
the STATCOM within the power system is performed 
in the MATLAB/Simulink environment using the 
Power System Blockset PSB. The performance of 
STATCOM scheme connected to the 230kV grid is 
evaluated. In addition, all details of the models 
implementation, the controls used, and the data for 
the test system are providing in the paper. 

Keywords: FACTS, STATCOM, Reactive Power 
Control. 

1. Introduction 

The STATCOM is the state of the art technology for 
dynamic control of reactive power. A STATCOM can 
be used for steady-state voltage control, to increase 
the stability of generator systems during transient grid 
fault and for dynamic reactive power generation. The 
STATCOM technology is based on a Voltage Source 
Converter connected by a filter and possibly a 
coupling transformer in shunt to the power system. A 
detailed review on literature about the STATCOM 
technology is given in [1]. Compared to the Static 
Var Compensator SVC, which consists of a thyristor 
controlled reactor, the STATCOM is often the 
preferred solution due to its dynamic reactive power 
generation capabilities, the smaller size and the 
ability to generate reactive power during low grid 
voltage condition. 

The commercial availability of gate turn-off thyristor 
(GTO) switching devices with high power handling 
capability, and the advancement of the other types of 
power semiconductor devices, such as IGBT’s have 
led to the development of fast controllable reactive 



power sources utilizing new electronic switching and 
converter technology. These switching technologies 
additionally offer considerable advantages over 
existing methods in terms of space reductions and fast 
effective damping [1]. Benefits provided by FACTS 
devices are twofold: they are capable of increasing 
power transfer over transmission lines and can make 
these power transfers fully controllable [2]. The last 
decade has witnessed a wide interest in design and 
application of FACTS devices, both from academic 
and industrial research. Several enhancement in 
semiconductor technologies and manufacturing have 
made power electronic devices capable of handling 
high power levels GTO, IGBT, IGCT commercially 
available [3]. Many devices have been proposed in 
the past within the FACTS family, some of them 
based on thyristor switched reactors like the Thyristor 
Controlled Series Compensator TCSC and the SVC 
[4], the Static Synchronous Series Compensator 
SSSC [5], the Unified Power Flow Controller UPFC 
[6] and the Superconducting. Magnetic Energy 
Storage systems SMES [7], employ power electronic 
converters. While others, like the STATCOM present 
in this paper. 

The need for more secure and efficient electricity 
supply systems management has given rise to 
innovative technologies in power generation and 
transmission. The combined cycle power station is a 
good example of a new development in power 
generation and FACTS as they generally known, are 
new devices that improve transmission systems [8- 
11]. In addition, the economical utilization of 
transmission system assets is of vital importance to 
enable electric utilities in industrialized countries to 
remain competitive and to survive. In developing 
countries, the optimized use of transmission systems 
investments is also very important to support industry, 
create employment and utilize efficiently scarce 
economic and industrial resources. FACTS is the new 
engineering technology that responds to these needs. 
It significantly alters the way power transmission 
systems are developed and controlled together with 
improvement in asset utilization, power system 
flexibility and effective system performance [12-17]. 
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Especially for the grid integration of wind farms the 
STATCOM technology has gained a lot of attention 
in the last years. The increased amount of renewable 
energy sources connected to the power system has 
forced transmission system operators to tighten their 
grid code requirements. A detailed overview on grid 
codes is given in [18]. A STATCOM for the grid 
integration of wind power plants in a weak loop 
power system can be found in [19]. In [20] and [21] 
the STATCOM is compared with a SVC to increase 
the stability of fixed speed induction generator wind 
turbines during grid voltage dips. By compensating 
the high amount of reactive power drawn by these 
generators after they have accelerated during a grid 
voltage dip can significantly increase the critical 
clearing time which is a criteria for the stability of the 
system. To perform an indirect torque control for the 
same kind of generators during grid voltage dips, the 
STATCOM is investigated in [22]. It can be shown 
that the torque of the machines can indirectly be 
controlled by the reactive power to reduce stress on 
the generator shaft to increase the gearbox lifetime. 

The coordinated control of a STATCOM and a wind 
farm with doubly-fed induction generators is shown 
in [23]. Most of these studies focus on the behaviour 
during three-phase grid voltage dips. This matches 
the requirements given by the grid codes, that often 
only consider balanced grid voltage conditions and 
only give short comments on the requirements during 
unbalanced grid voltage dips. But, the majority of 
grid faults are of unbalanced nature [24]. Such a 
negative sequence voltage can lead to undesired 
oscillations in current, STATCOM reactive and 
active power and DC link voltage and can make the 
STATCOM tripping, which should be avoided under 
all circumstances for power system stability reasons. 
Thus, safe control structures that can handle even 
significant voltage unbalances should be applied to 
the STATCOM. 

A STATCOM controller for handling unbalanced 
grid voltage with positive and negative sequence 
control loops can be found in [25]. Furthermore, the 
STATCOM currents can be injected into the power 
system aiming for different targets. The STATCOM 
can inject a purely positive sequence current or it can 
also inject a negative sequence current, which has 
been done in [26]. In [27] three control strategies 
under a given current rating of the converter are 
introduced. 

This paper proposes and validates models to 
accurately represent STATCOM regulates system 
voltage in transient by absorbing or generating 
reactive power of sample three-bus test system. The 
performance of STATCOM scheme connected to the 
230kV grid is evaluated. In addition, all details of the 
models implementation, the controls used, and the 
data for the test system are providing in the paper. 



2. STATCOM Model 

2.1 STATCOM Equivalent Circuit 

The STATCOM is a shunt device of the FACTS 
family using power electronics to control power flow 
and improve transient stability on power grids. The 
STATCOM regulates voltage at its terminal by 
controlling the amount of reactive power injected into 
or absorbed from the power system. The variation of 
reactive power is performed by means of a Voltage 
Sourced Converter VSC connected on the secondary 
side of a coupling transformer. The principle of 
operation of the STATCOM is explained on the Fig. 
(1). The active and reactive power transfer between a 
source V 7 and a source V 2 can be calculated by: 

VV 2 sin 8 

1 2 (i) 



p= 



X 



Where, V 7 represents the system voltage to be 
controlled and V 2 is the voltage generated by the VSC. 




Figure 1 . The active and reactive power transfer between sources 



In steady state operation, the voltage V 2 generated by 
the VSC is in phase with Vj ( 3=0 ), so that only 
reactive power is flowing ( P=0 ). If V 2 is lower than 
Vj, the reactive power Q is flowing from Vj to V 2 
(STATCOM is absorbing reactive power). While, if 
V 2 is higher than V 7 , the reactive power Q is flowing 
from V 2 to V 7 (STATCOM is generating reactive 
power). The amount of reactive power is given by: 



c V,(V t -Vi) 

X 



( 2 ) 



A capacitor connected on the dc side of the VSC acts 
as a dc voltage source. In steady state the voltage V 2 
has to be phase shifted slightly behind Vj in order to 
compensate for transformer and VSC losses and to 
keep the capacitor charged. Two VSC technologies 
can be used for the VSC: The first, VSC using GTO 
based square wave inverters and special 
interconnection transformers. In this type of VSC, the 
fundamental component of voltage V 2 is proportional 
to the voltage V^ c . Therefore, Vd c has to be varied for 
controlling the reactive power. The second, VSC 
using IGBT based PWM inverters. This type of 
inverter uses PWM technique to synthesize a 
sinusoidal waveform from a dc voltage source with a 
typical chopping frequency of a few kilohertz. This 
type of VSC uses a fixed dc voltage V^ c . Voltage V 2 is 
varied by changing the modulation index of the PWM 
modulator. 
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2.2 STATCOM Controls Circuit 

The control system consists of: A PLL which 
synchronizes on the three-phase primary voltage V h 
and the output of the PLL (angle 6= cot) is used to 
compute the direct-axis and quadrature-axis 
components of the ac three-phase voltage and 
currents labeled as V& V q or I d , I q on the Fig. (2). 
Measurement systems measuring the d and q 
components of ac voltage and currents to be 
controlled as well as the dc voltage V dc . An outer 
regulation loop consisting of an ac voltage regulator 
and a dc voltage regulator. Where: V is a positive 
sequence voltage, I is a reactive current, and X s is a 
slope or droop reactance. The output of the ac voltage 
regulator is the reference current I qre f for the current 
regulator. 

The output of the dc voltage regulator is the reference 
current I dre f for the current regulator. The current 
regulator controls the magnitude and phase of the 
voltage generated by the PWM converter ( V 2d , V 2q ) 
from the I dre f and I qre f reference currents produced 
respectively by the dc voltage regulator and the ac 
voltage regulator. The current regulator is assisted by 
a feed forward type regulator which predicts the V 2 
voltage output ( V 2d , V 2q ) from the Vj measurement 
(V Jd , V i q ) and the transformer leakage reactance. 




The STATCOM can be operated in two different 
modes: In voltage regulation mode, and other, in var 
control mode. When the STATCOM is operated in 
voltage regulation mode, it implements the following 
V/I characteristic as shown in Fig. (3). 




Figure 3. STATCOM V-I characteristic 



As long as the reactive current stays within the 
minimum and minimum current values (- I max , I ma x ) 
imposed by the converter rating, the voltage is 
regulated at the reference voltage V re f. However, a 
voltage droop is normally used (usually between 1% 
and 4% at maximum reactive power output). In the 
voltage regulation mode, the V/I characteristic is 
described by the following equation: 

V = V ref -X s I (3) 

Figure (4) shows the basic control structure of the 
STATCOM, based on a decoupled, d-q controlled 
strategy and the two levels, 24-pulse converter. The 
regulation slope k is defined as: 

hVcmax AF i max 

Icmax hmax 

Where AV cmax is the maximum line voltage drop at 
full capacitive and AV/ max is the maximum line over 
voltage at full inductive operation of the STATCOM. 
I cmax and Ii max are the corresponding STATCOM 
current ratings. 



Rcguinboa Slope 




V t 

Figure 4. Basic control structure of the STATCOM 



2.3 Dynamic Representation of the STATCOM 

A power system representation as seen by the 
STATCOM (including coupling transformer), and X eq 
is the equivalent Thevenin reactance of the power 
system as seen by the STATCOM, V eq is the 
equivalent Thevenin voltage and V, is the STATCOM 
bus voltage. For the STATCOM capacitive operation 
we have: 

Vt = V eq + I q x eq ( 5 ) 

Figure (5) shows the functional block diagram 
representation of the STATCOM when operating in 
capacitive mode. The STATCOM can be modeled as 
a reactive current source with transportation lag, T q . 
The voltage regulator and the feedback voltage 
measurement are also modeled with transportation 
lags of T l9 and T 2 , respectively. 




Figure 5. Block diagram representation of the STATCOM 
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2.4 STATCOM Mathematical Equations 

The differential equations for a STATCOM 
summarized are: 



d L 

dt 
d i ' 



dt 

di c S 

dt 

dv 



' = A 0 dJ^i R Ls-( V aS- V as)) 



DC 



dt 



3l 

b 



v DC 



+ 1 



DC 



( 6 ) 

(7) 

( 8 ) 
(9) 



Where: i DC and v‘ as , v bs , v‘ cs are the dc side current 
and converter output phase voltages respectively. 
Now, (i as + i bs + i cs ) = 0 and if we assume 

(v as + V bs + V rl ) = 0 then it implies that 
(v' +vi + v‘ )=0. 

\ as bs cs ) 



Then we can write the converter output voltages as: 



Ks = V DC S 6 W 


(10) 


V l = V DC S 6 W 


(11) 


V ‘cs = V DC S e(0 


(12) 


Where S£ ,S b 6 , S c 6 are 


the switching function for a 



six-pulse STATCOM, In the steady state: 






t + - 



V 



3 j 



= s: 



t + ■ 



An 

j 



(13) 



If we assume the converter to be lossless, 

(14) 



V DC l DC = V as l as+ V bs l bs+ V cs l cs 



Therefore, by substitute equation (11) to (13) into 
equation (14) 



*dc — i as S 6 + hs$ 6 + i a S 6 



(15) 



Cleary, the dynamical equation are time-variant and 
difficult to analysis for a varying firing angle .If the 

switching function ,S£ , S£ are approximated by 

their fundamental since component (neglecting 
harmonics) we get: 

2 

V as = v DC — sin (co 0 t + a + 0) (16) 

n 

2 

v l bs = v DC — sin(^ 0 + a + 6 - 2^/3) (17) 

n 



v = v 



DC 



71 



-sin(&> 0 + a + 0 + 2x/3) 



( i as 5m(o) 0 + a + O') 

+i bs sin((jL> o + a + 6 — 2n/3 ) 
+i as sin(cu 0 + oc + 9 + 27r/3), 



(18) 



(19) 



The approximation made above is valid if the 
harmonic component are small, compared to the 
fundamental. For a six pulse STATCOM harmonic 
components are substantial. However, if higher pulse 
numbers are used, harmonic can be reduced. The 
above analysis for a six pulse STATCOM can be 
easily extended to higher pulse number where the 
approximations are justified. From the point of view 
of control system design, we can neglect the 
harmonics though detailed simulation is necessary to 
validate any results. The STATCOM bus (phase) 
voltages are given by: 

„ = ^|2pV, sin(ffl 0 t + 6) (20) 

v»,=VV3V, sin(ty 0 f -\-d-2nj2) (21) 

sin(tfy + 0 + 2^/3) (22) 

Where V s is the line-to-line rms voltage and 0 is the 
angle of the STATCOM bus. Note that a is the angle 
by which the fundamental component of converter 
output voltages leads the STATCOM bus voltages. 



We use Kron's transformation to the D-Q-0 variables: 



v bs 



= [C.] 



V DS 



QS 



(23) 



i 



L os 



Where: 




cos(^) 
cos(/y 0 f — 2 tt/ 3) 
cos(tfy + 2 tt/3) 



sin(fty) 
sin(co 0 t -2^/3) 
sin(#y + 2/r/3) 



l/V2 

l/V2 

l/V2 



The ac phase voltages can also be transformed in a 
similar manner. Since [c k ] ' = \c k J the transformation 

is power invariant, we can ignore i 0 because it is 
always zero. The equation in the D-Q variables are 
time invariant and are given by: 

di r 



l ds 



dt 

di, 



QS 



dt 



di 



DC 



dt 





l DS 


r “| 




l QS 


+ ^\: c 
L J 




Jos _ 





(25) 



Where: 



[*i]= 



- Roj () 



-a>n 



co 0 ksin(a+d) 



cd^ 



- Roj () co 0 kcoi(a+d) 



co 0 k sin(a + 6) co 0 k co i(a + 6) 



b c R p 
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[K 2 ] 



~ a) o 

0 

0 



0 

~ G) o 

0 



6 = tan 1 ( v DC /v QS ) angle of STATCOM bus voltage, 

v s = ^{ v ds + v qs ) Ma g nitude of STATCOM bus 
voltage, OC is the angle by which the converter 
output voltage leads the bus voltage, and k = Vb/^ 



2.5 Decoupled Current Control System 

The decoupled control system is based on a full d-q 
decoupled current control strategy using both direct 
and quadrature current components of the 
STATCOM ac current. The decoupled control system 
is implemented as shown in Fig. (6). A PLL which 
synchronizes on the positive sequence component of 
the three phase terminal voltage. The output of the 
PLL is the angle 0 which used to measure the direct 
axis and quadrature axis component of the ac three 
phase voltages and current. The outer regulation loop 
comprising the ac voltage regulator which provide the 
reference current I qre f for the current regulator which 
is always in quadrature with the terminal voltage to 
control the reactive power. The voltage regulator is a 
proportional plus integral PI controller with K p and K,. 
The current regulator is also PI controller with K p and 
Kj. The PLL system generates the basic 
synchronizing signal which is the phase angle of the 
transmission system voltage V s , 6 and the selected 
regulation-slope, k , determines the compensation 
behavior of the STATCOM device. 




Figure 6. STATCOM d-q decoupled current control system 



3. Studied System 

The studied system is sample three-bus study system 
with the STATCOM located at bus B 2 consists of 
three phase ac source has supple one transmission 
line by 230kV and these transmission line supple 
three loads. The first load is inductive load (l+j0.8)pu, 
the second load is inductive load (0.7+j0.5)pu and the 
third load is change between two value the first value 



capacitive load (0.6-j0.4)pu and the second value is a 
dynamic load as shown in Fig. (7). 




Figure 7 Three-bus study system with the STATCOM located at 
Bus B 2 



The three-phase source implements a balanced three- 
phase voltage source with internal R-L impedance, 
and a transmission line parameter model with lumped 
losses are based on the Bergeron's traveling wave 
method used by the Electromagnetic Transient 
Program (EMTP). The lossless distributed LC line is 
characterized by two values, the surge impedance 
Z c WlC an d the phase velocity v=l/VLC The 

model uses the fact that the quantity e+Zi (where e is 
line voltage and i is line current) entering one end of 
the line must arrive unchanged at the other end after a 
transport delay of T — d / V , where d is the line 
length. By lumping R/4 at both ends of the line and 
R/2 in the middle and using the current injection 
method, the following two-port model is derived as 
shown in Fig. (8). 



j* it 




Where: 

Z=Z C (R/A) 

,, Z C -{R/ 4) 

Z C +(R/ 4) 

Z c =jL/C 

T=d^[LC 

The three-phase transformer implements a three- 
phase transformer using three single-phase 
transformers. The loads implement a three-phase 
balanced load as a parallel combination of RLC 
elements. At the specified frequency, the load 
exhibits constant impedance. The three-phase 
dynamic load implements a three-phase, three-wire 
dynamic load whose active power P and reactive 
power Q vary as function of positive-sequence 
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voltage. The load impedance is kept constant if the 
terminal voltage V of the load is lower than a 
specified value V min . When the terminal voltage is 
greater than the V min value, the active power P and 
reactive power Q of the load vary as follows: 

/ T7 \ n P A rr, \ 



p (s)= p , 77- 



Q( s )=Qc ir 




yx -r ^ pl o j 

(l + T p2 s ) 


(26) 


(i+v) 

(i+v) 


(27) 



Where: V 0 is the initial positive sequence voltage, P Q 
and Q 0 are the initial active and reactive powers at the 
initial voltage, V is the positive-sequence voltage, n p 
and n q are exponents (usually between 1 and 3) 
controlling the nature of the load, T pl and T p2 are time 
constants controlling the dynamics of the active 
power P, and T ql and T q2 are time constants 
controlling the dynamics of the reactive power Q. 
Table 1. Shows the three bus studied system with 
STATCOM parameters. Also, phase shift pattern for 
24-pulse converter as in Appendix. 



Table 1. Three bus studied system consist of two loads with 
STATCOM parameters. 



Three phase ac source 

Rated Voltage: 1.03p.u 
Short circuit capacity: 10000 
MV A 

Frequency: 60Hz 
Base Voltage: 230 kV 
X/R ratio: 8 


STATCOM: 

Primary Voltage: 138 kV 
Secondary Voltage: 15 kV 
Nominal Power: 100MVAR. 
Frequency: 60 Hz 
Equivalent Capacitance: 

750pF 


Power Transformer(Y/A): 

Primary Voltage: 230 kV. 
Secondary Voltage: 33 kV. 
Frequency: 60 Hz 
Rated Power: 300 MVA 
Magnetization Resistance: 
500 p.u., and Magnetization 
Reactive: 500 p.u 


Coupling Transformer 

(Y/Y): 

Primary Voltage: 138 kV 
Secondary Voltage: 230 kV 
Frequency: 60 Hz 
Rated Power: 100 MVA 


Capacitor Bank (dc): 

Total Capacitance: 10 mf 
Rated dc Voltage: 2 kV 


Transmission Line 

Resistance: R T l- 0.05 p.u. 
Reactance: X T c 0.2 p.u. 


GTO Switches 

Snubber Resistance: 100 KX2 
Snubber capacitance: inf. 
Internal Resistance: 0.1 Q 
No. of Bridge arm: 3 


Three Phase Loads 
Load 1 (1.0+j0.8) p.u 
Load 2 (0.7+j0.5) p.u 
Load 3 (0.6-j0.4) p.u 



4. Results and Discussions 



The system study radial power system is subjected to 
load switching at Bus B 3 . At starting, the source 
voltage is such that the STATCOM is active. It does 
absorb provide reactive power to the network. The 
capacitor bank is recharged to 1.0 pu voltage. The 
network voltage, V g , is 1.02 pu and only capacitive 
load 3 with (P=0.6pu and Q-OA pu) (at rated voltage) 
is connected at load Bus B 3 and the STATCOM B 2 
bus voltage is 1.2pu for the uncompensated system 
and the transmitted real and reactive power are 
P L =0.7pu and Q L =-0Apu. 



The results can summarize for the four cases are the 
following load excursion sequence: 



Case 1: At t= 0.5 sec load 1 inductive load ( l+j0.8 pu) 
injection, at t-1 sec load 2 inductive load (0.7+j0.5) 
injection, at t=1.5 sec load 3 capacitive load and load 
2 inductive load rejection. 

Case 2:At t=0.5sec load 2 inductive load ( 0.7+j0.5pu ) 
injection, at t=l sec load 1 inductive load (l+j0.8) 
injection, at t=1.5 sec load 3 capacitive load and load 
2 inductive load rejection. 

Case 3: At t=0.5 sec load 1 inductive load ( l+j0.8 pu) 
injection, at t-1 sec load 2 inductive load (0.7+j0.5) 
injection, at t-1. 5 sec load 3 Dynamic load and load 
2 inductive load rejection. 

Case 4:At t=0.5sec load 2 inductive load ( 0.7+j0.5pu ) 
injection, at t=l sec load 1 inductive load (l+j0.8) 
injection, at t=1.5 sec load 3 Dynamic load and load 
2 inductive load rejection. 

Figure 9(a) is represented case 1, the terminal voltage 
without STATCOM is 1.2 pu starting and 0.96 pu 
after injection inductive load 1, 0.85 pu after injection 
inductive load 2, and near unity after rejection of both 
capacitive (load 3) and inductive (load 2) loads. 
Whenever, with STATCOM the sequences of the 
terminal voltage are 1.02 pu, 1.0 pu, 0.97 pu, and 
unity respectively. The active and reactive powers of 
the transmission line are starting (0.9-j 1 .3) pu without 
STATCOM and are (1.43+jO.l) pu after inductive 
load 1 injected, (1.5+j0.75) pu after inductive load 2 
injected, and back to (0.7-j0.25) pu after rejected of 
both inductive load 2 and capacitive load. While with 
STATCOM the sequence of the active and reactive 
powers on the transmission line are: (0.7-j0.4) pu, 
(1.5-j0.1) pu, (2.0+j0.2) pu, and (0.7-j0.25) pu 
respectively as shown in Fig. 9(e, f). A DC voltage 
and currents of the VSC STATCOM as shown in Fig. 
9(b, c, d) respectively. 

From Fig. (10) to Fig. (12) are depicts the terminal 
voltage, DC voltage, voltage and current of the VSC 
STATCOM, active and reactive power in the 
transmission line, all these in three cases (case 2 to 
case 4) respectively. Figure (13) discuses the reactive 
power supply and absorbing by STATCOM for 
different cases. When load 1, and load 2 are inductive 
load and load 3 is capacitive load as shown in Fig. 
10(a). While the reactive power supply and absorbing 
by STATCOM when load 1, and 2 are inductive load 
and load 3 is Dynamic load as shown in Fig. 13(b). 
STATCOM absorbed reactive power 0.8pu at starting 
and delivered 0.2 pu after injection inductive load 1, 
0.82pu after injection inductive load 2, and near zero 
after rejection of both capacitive load 3 and inductive 
load 2 loads. 
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(c) Voltage of the (VSC) STATCOM. 



1.2OE+O0 

§ iooe+oo ■ 

5 SjOOE-OI ■ 

g 4.00C-01 
s Z.00E-01 - 

^ n nftc - HA 




. . . . 






-with autSTATCOM 

— with STATCOM 


' 


Jl . 




1 " " 






























H 




























_ . _ < 




_ „ „ < 










. . . . 










1 




1 








" - " 1 










r 












V 






UUIPL+W 
“ 1 


> 0 


2 o 


4 O 


6 O 


A i 


1 


2 1 


J4 1 


,6 1 


a : 


Time4sec} 





(d) Current of the (VSC) STATCOM. 
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(f) Reactive power in the transmission line. 

Figure 9. Transient performance of the STATCOM when load 1, 2 
inductive load and load 3 capacitive load (case 1). 
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(b) DC voltage of the STATCOM 
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(f) Reactive power in the transmission line. 

Figure 10. Transient performance of the STATCOM when load 1, 
2 inductive load and load 3 capacitive load (case 2). 
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(a)Terminal voltage at bus 2 




(b) DC voltage of the STATCOM 




(c) Voltage of the (VSC) STATCOM. 
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(d) Current of the (VSC) STATCOM. 
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(e) Active power in the transmission line. 
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(f) Reactive power in the transmission line. 

Figure 11. Transient performance of the STATCOM when load 1, 
2 inductive load and load 3 Dynamic load (case 3). 
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(b) DC voltage of the STATCOM 
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(d) Current of the (VSC) STATCOM. 
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(e) Active power in the transmission line. 
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Figure 12. Transient performance of the STATCOM when load 1, 
2 inductive load and load 3 Dynamic load (case 4). 
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Figure 13. Reactive power supply and absorbing by STATCOM 
for different cases 



5. Conclusions 

Numerical simulations are reported for a STATCOM 
based on the proposed converter, showing simplicity 
both in control and converter topology, as well as 
erective tracking performance of the sliding controller. 
GTO voltage source converter of STATCOM device 
model is presented. These full descriptive digital 
models are validated for voltage stabilization reactive 
compensation in a sample three-bus study system. 
The performance of the STATCOM when injection 
and rejection loads. The maximum load ability and 
voltage profiles of power system can be significantly 
improved by STATCOM. When hitting the current 
limit, the related system buses will lose the 
controllability and bus voltages decreased rapidly 
along with the loading level. 
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Appendix 

A phase shift pattern for 24-pulse converter are: 

PST +7.5° Necessary to eliminate the 24-pulse harmonics 
+3.75° Necessary to eliminate the 48-pulse harmonics 

Total +11.25° Winding turn rate 1: tan (11.25°) 

Driver -7.5° Necessary to eliminate the 24-pulse harmonics 
-3.75° Necessary to eliminate the 48-pulse harmonics 



v aM 2 ( f )i = 2[V m sin( wt + 30°) +V abn sin(l Iwt + 195°) 
+ V abl3 sin(l 3wt + 255° ) + . . . 

+ V ^23 sin(23 wt 4- 60° ) 

+ V a ^25 sin(25wf + 120+ + • • •] 

PST -7.5° Necessary to eliminate the 24-pulse harmonics 

+3.75° Necessary to eliminate the 48-pulse harmonic 



Total -3.75° Winding turn rate 1: tan (3.75°) 

Driver +7.5° Necessary to eliminate the 24-pulse harmonics 

-3.75° Necessary to eliminate the 48-pulse harmonics 



Total +3.75° 

The resultant output voltage generated by the second 12-pulse 
converter is: 

V abl2 (t) 2 = 2[ V ab] sin( wt + 30" ) + V ab] , sin(l Iwt + 15°) 

+ sin(13wT + 75 ) + 

+ V ab23 sin(23wf + 60°) 

+ V ab25 sin(25wf + 120°) + ---] 
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Total -11.25° 

The resultant output voltage generated by the first 12- pulse 
converter is: 
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Resume 

Un systeme dynamique hybride est decrit par un 
ensemble de variables continues et un ensemble 
d’ evenements discrets interagissant mutuellement. La 
realite impose en outre de prendre en compte les 
defaillances des composants ou les incertitudes sur la 
connaissance du systeme. L’accroissement incessant 
de la complexity des systemes automatiques a fait en 
sorte que la tache de la surveillance devienne de plus 
en plus moins aisee. De ce fait, les systemes peuvent 
fonctionner en plusieurs modes. Certains de ces 
modes correspondent a un fonctionnement normal et 
d’autres represented des modes defaillants. Dans cet 
article, nous nous interessons a 1’ evaluation de la 
probabilite d’ occurrence des evenements defaillants 
dans un systeme dynamique hybride. Nous proposons 
pour cela, une methode d’ evaluation basee sur la 
connaissance a priori du systeme, au moins du point 
de vue de l’etat de ses composants en fonctionnement 
normal (ferme ou ouvert, eteint ou allume,... etc) 
Afin de valider notre proposition, nous l’avons 
appliquee au benchmark utilise par differents auteurs 
pour le calcul de la fiabilite dynamique [1] . Notre 
objectif est done d’evaluer la probabilite d’ occurrence 
des deux evenements redoutes du cas test : 
l’assechement et le debordement d’un reservoir. Cet 
exemple a pour principal interet de mettre en lumiere 
les problematiques de la surete de fonctionnement des 
systemes. 

Abstract 

A dynamic hybrid system is described by a set of 
continuous variables and a set of discreet events 
interacting mutually. The reality imposes to take into 
account the failures of constituents (components) or 
the uncertainties on the knowledge of the system. The 
ceaseless increase of the complexity of the automatic 
systems made so that the task of the surveillance 
becomes more and more less easy. Therefore, the 
systems can work in several modes. Some of these 
modes correspond to a normal functioning and the 



others represent failing modes. In this article, we are 
interested in the evaluation of the probability of 
occurrence of the failing events in a hybrid dynamic 
system. We propose for it, a valuation method based 
on the knowledge in priori of the system, at least from 
the point of view of the state of its components in 
normal functioning (closed or opened, put out or lit, 
etc.). To validate our proposition, we applied her to 
the benchmark used by various authors for the 
calculation of the dynamic reliability. Our objective 
thus is to estimate the probability of occurrence of 
both events dreaded by the test case: the drying out 
and the overflowing of a reservoir. 

This example has for main interest to bring to light 
the problems of the safety of functioning of the 
systems. 

Mots-Cles : Diagnostic, systemes dynamiques 

hy brides, surveillance, surete de fonctionnement. 

Keywords: Diagnosis, dynamic hybrid systems. 

Reliability 

1. Introduction 

NOUS representons un systeme, d’un point de vue 
diagnostic comme etant un cycle de vie qui est 
interrompu par 1’ occurrence d’un evenement qui fait 
basculer le mode de fonctionnement du systeme d’un 
etat normal a un etat de fonctionnement defaillant. En 
adoptant cette representation, le concepteur peut 
aisement remarquer que son apport ne peut se 
manifester qu’au niveau du bloc detection pour 
ensuite reconfigurer le systeme. L’ occurrence de 
T evenement, on ne peut plus aleatoire. L ‘objectif de 
cette etude est precisement de denombrer et 
caracteriser les etats defaillants ou ce que nous 
appelons etats redoutes et d’evaluer la probabilite 
d’ occurrence des evenements qui menent vers l’etat 
defaillant. 

Un systeme est un ensemble de composants et le 
composant, dans notre cas, est considere comme une 
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entite represente par un nombre binaire (bit). Le 
fonctionnement du systeme est divise en etapes 
temporelles. Chaque etape est consideree comme telle 
avec tous ses parametres. On considere une etape de 
fonctionnement comme un sous- etat de 

fonctionnement qui peut etre un etat normal comme il 
peut etre un etat defaillant. Dans le cas le plus simple, 
le systeme peut etre considere comme une seule etape. 
Les parametres de l’etat sont les valeurs des bits 
representant les composants du systeme (1 ou 0 par 
exemple ferme, ouvert ; eteint, allume ...), on 
n’entend pas ici defaillant ou non. Le sous-systeme 
avec 1’ etape temporelle correspondante en question 
est alors une combinaison des bits representatifs des 
composants. Une etape faisant intervenir N 
composants aura 2 N combinaisons representant 2 N 
sous-etats possibles (defaillants et normaux). Aussi, 
chaque combinaison nous informe de l’etat de chaque 
composant et indirectement, elle nous informe aussi 
sur les composants qui sont a l’origine de l’etat 
defaillant. Une fois les sous-etats sont simules, les 
resultats sont analyses et les modes defaillants sont 
reperes. Apres quoi, un calcul de probabilite est 
effectue. 

4- Decomposer le fonctionnement total du systeme 
en etapes de fonctionnement ou les composants 
intervenant ont un seul etat. 

4- Reperer les composants intervenant dans chaque 
etape 

4- Representer chaque composant par un bit B {1,0} 
1 : Etat demande et 0 : Etat contraire 
4- Representer V etape par un nombre binaire dont le 
nombre de bits est egal au nombre de composants 
intervenant dans 1’ etape en question 
4- Exprimer les combinaisons representatives de 
1’ etape (sous-etats de fonctionnement relatif a 
1’ etape) 

4- Analyser les resultats 

4 Calculer la probabilite d’ occurrence de chaque 
etat defaillant 

2. Le cas test 

La representation d’un systeme dynamique d’un point 
de vue diagnostic est schematisee sur la figure. 1. II est 
considere comme etant un cycle de vie qui est 
interrompu par V occurrence d’un evenement qui fait 
basculer le mode de fonctionnement du systeme d’un 
etat normal a un etat de fonctionnement defaillant Le 
cas test consiste en un reservoir contenant un liquide 
dont le niveau h doit etre maintenu a l’aide d’une 
pompe principale PI, d’une pompe de secours et 
d’une vanne de vidange V. 

Chacun des trois composants est commande par une 
boucle de controle contenant un detecteur de niveau 
(figure 2). 




Figure. 1. Systeme de regulation de niveau 



Signal de niveau Signal de niveau 




Figure. 2. Etat global d’un systeme dynamique 



La mission a remplir est de maintenir le niveau du 
liquide constant dans un intervalle ([6,8]). Si le niveau 
est inferieur a 6, la commande ferme la vanne et 
ouvre la pompe de secours et si le niveau est superieur 
a 8, les pompes sont arretees et la vanne est ouverte. 
Deux situations extremes peuvent se produire : 
l’assechement et le debordement. Ces deux cas 
arrivent quand la commande ne peut plus agir sur les 
composants du systeme qui deviennent alors 
defaillants (1 ou plusieurs composants). On dit alors 
que le systeme est dans un etat defaillant voir 
dangereux. Les trois composants la pompe PI, la 
pompe P2 et la vanne V sont mutuellement 
independants et non reparables. 

La variable continue pour le systeme est le niveau du 
liquide h lequel est fonction de l’etat des composants. 
Ainsi, 1’ equation differentielle pour le systeme est 
donnee par : 
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^ = K«) w 

OU U = ( U U p2, U y ) 



Pour le reservoir nous avons les etats ; 

- N_n : niveau normal du reservoir ( 6 ^ h 8) 

- N_ass : niveau d’assechement ( h ^ 4) et 

- N_deb : niveau de debordement ( h ^ 10 ). 



Avec V c 



{? 



si c est OFF 
si c est ON 



avec y (u ) = (u pl + u P2 _u v )D 



( 2 ) 



ou D est le debit des elements PI, P2 et vanne. 

L’ equation generalisee (1) reflete les differents modes 
operatoires possibles du processus. Elle fait apparaitre 
1’ influence des phenomenes discrets sur V evolution 
du processus au travers des termes u pi, u p2, u v- Ces 
derniers peuvent prendre, dans le cas de cet exemple, 
la valeur 1 si le composant est actif ou bien si une 
defaillance le bloque dans la position ouverte ou 
active, et la valeur 0 dans le cas contraire, comme 
1’ exprime V equation (2). En conditions nominates, le 
debit de PI est egal au debit de P2 et de V. Ainsi D = 
1,5 m 3 h'l pour PI, P2 et V. 

Au temps t = 0, le niveau du liquide h = 7 m, la 
pompe PI fonctionne, P2 est a V arret et la vanne V est 
ouverte. 

Les lois de commande du systeme qui definissent 
l’etat des pompes et de la vanne en fonction du niveau 
de liquide sont donnees dans le tableau 1 suivant : 



Tableau 1: Les lois de commande du systeme 



Niveau 

h(m) 


Pompe PI 


Pompe P2 


Pompe P3 


h<6 


Ouverte 


Ouverte 


Lermee 


6<h>8 


Ouverte 


Lermee 


Ouverte 


h>8 


Lermee 


Lermee 


Ouverte 



Le systeme ainsi defini peut etre considere comme un 
automate hybride qui prend en compte les differents 
modes continus de fonctionnement du systeme et le 
passage de l’un a V autre sur V occurrence des 
evenements determinates qui sont produits par 
franchissement de seuils des variables continues. Les 
automates elementaires sont donnes sur la figure 3. 

Les evenements associes a PI, P2, P3 sont: 

- a_Pl et a_P2 : arret des pompes PI et P2 

- d_Pl et d_P2 : demarrage des pompes PI et P2 

- o_V : ouverture de la vanne et 

- f_V : fermeture de la vanne 

Nous distinguons pour les composants PI, P2 et V les 
etats suivants : 

- ON-PI, ON-P2 et ON-V : pompes actives et 
vanne ouverte, 

- OLL-P1, OLL-P2 et OLL-Y : pompes en arret et 
vanne fermee, 



L’ automate representant la commande est donne sur 
la figure 4. Les lois de controle de la commande: 

- initial : CO - PI active, P2 arretee et V ouverte, 

- si h < 6 : C2 -► C3, PI active ; C3 C4, P2 
active et C4 -► C5, V fermee, 

- si h > 8 : C6 - C7, PI arretee ; C7 - C8, P2 
arretee et C8 — * Cl, V ouverte, 




(a) Automate de la pompe PI 




(b) Automate de la pompe P2 




d_V 



a_V 
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(d) Automate du reservoir 



La simulation d’une panne a la 9 eme ’iteration (vanne 
bloquee en fermeture) le niveau du liquide augmente 
indefiniment : debordement et l’indicateur de cet etat 
defaillant se met a 0 (figure 6 ). La simulation d’une 
autre panne (la pompe PI n’ arrive pas a s’ouvrir et la 
pompe P2 est bloquee en fermeture, la vanne etant 
ouverte), le niveau descend au plus bas niveau et 
l’indicateur de l’etat assechement est passe a 0 . II 
faut noter que les indicateurs des deux etats 
assechement et debordement ont ete initialises a 1 . 



Figure. 3. Automates a etats finis des trois composants 




3. Application de la methode 

Detection des etats defaillants 

Le fonctionnement total du systeme peut-etre decrit 
par une seule etape ou bien un seul etat a savoir 
maintenir le niveau du liquide h dans l’intervalle [ 6 , 8 ]. 
La simulation du systeme sans occurrence 
d’evenements intempestifs donne la courbe suivante 
sachant que le niveau initial est h=7m. Le systeme vu 
comme 9 a, maintient le niveau du liquide constant 
(figure 5). 



Volume h du reservoir Bat debordement Bat assechement 



8 

7.5 
7 

6.5 
6 




Figure. 5. Fonctionnement nominal du systeme 
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a) Assechement b)Debordement 

Figure. 6 . Etats defaillants du systeme 



Calcul de Probability 

La figure 7 represente les differentes combinaisons 
possibles (sous-etats). Le systeme demarre avec l’etat 
initial normal (101 : Pompe PI ouverte, Pompe P2 
fermee et la vanne V ouverte). A t=9, il ya 
V occurrence de Vetat defaillant (un , deux ou les 
trois composants ne sont pas dans Vetat attendu). 
Pour eviter toute ambiguite avec le role de la 
commande, celle-ci est inhibee a 1 ’ occurrence de 
l’evenement redoute. Par ailleurs, on ne s’interesse 
pas ici a la cause physique de la defaillance. Celle-ci 
peut etre due, en effet, a une action externe (fermeture 
ou ouverture intempestive), au blocage du composant 
dans son etat (ouverture ou fermeture) ou tout 
simplement cette defaillance est la consequence de la 
duree de vie du composant (usure). A partir de l’etat 
normal, on s’interesse alors a toutes les possibility de 
defaillances. Ces modes de fonctionnement sont 
resumes dans ce qui suit, a noter que seul l’etat 101 
est l’etat de fonctionnement nominal. 

1/ 000: La pompe PI est fermee, la pompe P2 est 
fermee et la vanne est fermee. Etat stationnaire dans 
l’etat initial. 
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2/ 001 La pompe PI fermee, la pompe P2 fermee et 
la vanne ouverte Assechement 
3/ 010 la pompe PI est fermee, la pompe P2 est 
ouverte et la vanne est fermee Debordement 
4/ 011 la pompe PI est fermee et la pompe P2 est 
ouverte, la vanne est ouverte. La pompe de secours 
prend le relais de la pompe Pl. Etat normal (de 
secours) 

5/ 100 la pompe PI est ouverte la pompe P2 est ferine 

et la vanne fermee Debordement 

7/ 110 la pompe PI ouverte la pompe P2 ouverte et la 

vanne fermee Debordement 

8/ 111 la pompe P 1 ouverte la pompe P2 ouverte et la 

vanne ouverte Debordement 



1000 

0 

-1000 

I 

8 




001 



Etat Normal Assechement Debord 

ement 

Probability 0,375 0,125 0,5 

Tableau 3 : Valeurs de probability de chaque etat 

Les resultats obtenus montrent que la probability de 
l’etat assechement est de 0,125 et la probability 
d’ occurrence de l’etat debordement est de 0,5. Si nous 
observons les valeurs de probability obtenues dans les 
etudes relatives a la fiabilite dynamique a t=1000 
Heures [1] [2], les resultats sont tres proches et si on 
augmentait le temps a plus de lOOOh, les valeurs 
d’ occurrence des deux evenements s’approcheraient 
de plus en plus des valeurs que nous avons trouve 
jusqu’a se stabiliser a celles-ci. 

Tableau 4 : Probabilites d’acces aux etats redoutes 



2000 



7 
6 

l 

8 
7 
6 

l 





Temps(h) 


Debordement Assechement 




Oil 




1000 


PMDPM RdP ASH PMDPM RdP ASH 

0,486 0,486 0,475 0,118 0,118 0,117 





101 




0 10 20 30 0 10 20 30 

Figure. 7. Niveau du liquide h pour tous les modes 



4. Conclusion 

Dans ce travail, nous avons propose une strategic 
pour le diagnostic des systemes dynamiques hy brides. 
Nous avons, ainsi, introduit une approche conduisant 
a fournir une aide pour V identification et la 
localisation des modes defaillants ainsi que le calcul 
de probability de leurs occurrences. Notre approche 
nous a permis d’evaluer la probability d’ occurrence 
des evenements redoutes du cas test. Les resultats 
trouves sont tres probants. La simulation a ete faite 
sur matlab. La prochaine etape est d’appliquer cette 
methode sur des exemples encore plus complexes. 



Les resultats de la simulation ont monte que tous les 
modes defaillants ont ete detectes. Ainsi, le mode 
debordement apparait dans 4 cas, le mode 
assechement 1 fois et l’etat normal 3 fois (Tableau 2). 



Etat normal 


Assechement 


Debordement 


000 


Oil 


010 


Oil 




100 


101 




110 






111 



Tableau 2 : Tableau recapitulatif des etats 



Le mode correspondant a la combinaison 000 ou tous 
les composants sont fermes, le niveau du liquide ne 
varie pas de l’etat initial. Pour le mode (Oil), la 
pompe PI est fermee et la pompe P2 prend le relais, 
cet etat correspond a l’etat normal et le niveau du 
liquide est constant. 

Le calcul de probability d’ occurrence des modes 
defaillants donne explicitement les valeurs (Tableau3) 
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Abstract 

This paper suggests a strategy for analysis and design of 
a digital robust control system, capable of suppressing 
the effects of the variations of its parameters. The 
methodology used for the robust controller design is 
compensation with two degrees of freedom that brings as 
a result a desired system’s performance and stability. The 
control system is analyzed in terms of stability with the 
aid of the D-partitioning method before and after being 
compensated. The system’s performance is also 
examined for robustness in the time domain. The analysis 
shows a considerably improved performance in terms of 
system’s robustness. The suggested robust controller 
build by two microcontrollers can be applied to any 
complex digital control system with variable parameters. 

Keywords: Variable parameters, D-partitioning, 

Discrete-time, Robust controller, Microcontrollers; 

1. Introduction 

The objective of the paper is to propose a strategy for 
analysis and design of an optimal robust controller that 
can be applied to digital control systems with variable 
parameters. It is common knowledge that the problem of 
digital compensators of the type of PI or PID controllers 
using difference equations is already resolved in 
literature, whereas this research is focusing mainly on the 
challenges of the digital robust control design. The 
technique employed for the robust controller design is 
achieved by applying series compensation with a unity 
feedback, combined with a forward compensation. The 
controller, operating in two degrees of freedom, enforces 
a desired system damping, margin of stability and time 
response. The designed robust controller suppresses the 
effects of parameter variations so that any parameter 
uncertainties would affect insignificantly the final system 
performance. The robust analysis and control of the 
digital system are discussed in the discrete-time domain. 
The lack of a comprehensive, simple and user friendly 
analysis tool for digital control systems and further the 
absence of a robust controller that can make a digital 
system insensitive for any parameter variation within 
specific limits are the motivating factors for this research. 
Following some initial ideas of Neimark’ s D-partitioning 



[1, 2], the method is further advanced in this research in 
terms of its application for digital control systems. The 
employment of the D-partitioning enables a quick and 
convenient determination of the regions of stability in 
case of variation of system’ s parameters and its effect on 
the system’s stability. The D-partitioning method is 
applied for analysis of digital system’s stability before 
and after implementing the designed robust controller. 
The employment of MATLAB enables automatic plotting 
and analysis of the regions of stability. Further analysis in 
terms of the system’s time response in the discrete-time 
domain is showing the efficiency of the designed digital 
robust controller. 

In this research, a simulation model, representing a 
complex real-time plant of high order is preliminary 
reduced and approximated to a third order system 
following a well known approximation criterion [3]. 
Although the complexity of the examined simulation 
model is significantly lower than the real-world plant, the 
achieved results can be implemented to any advanced 
system. In reality, most control systems contain elements 
that drift or vary with time and the initial correct 
determination of the exact parameter values is often quite 
challenging and problematic [3,4]. 

The shortcomings of a number of known optimization 
methods are the requirements of very accurate models of 
the systems. Optimizing systems with wrongly 
determined models often makes their performance worse 
rather than better. Alternatively, the suggested in this 
research analysis and design methods are to a certain 
extent insensitive to parameter drift. The applied 
compensation technique with two degrees of freedom, 
used for the design of the digital robust controller, even 
allows its operation in cases of simultaneous variation of 
a number of uncertain parameters, bringing the system to 
the desirable performance. 

2. The Concept of Digital Robust Control 

Systems 

A usually known block diagram of digital control system, 
shown in Figure 1, illustrates the interaction between its 
components. The interface at the input of the digital 
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computer is an analogue-to-digital converter (ADC). It 
converts the continuous-time error signal E(s) into a 
digital binary format signal that can be processed by the 
computer. At the computer output a digital-to-analogue 
converter (DAC) is required to converter the computer 
binary signals into continuous-time signal format that is 
necessary to drive the plant. A clock synchronizes the 
operation of the ADC, digital computer and DAC. The 
digital computer in combination with its interface, ADC 
and DAC, operates as a digital compensator within the 
system. A sensor, included in the negative feedback of 
the system, is representing a transducer that converts the 
plant’s variable under control into a continuous signal of 
a desired format [5, 6]. In the following discussions on 
the system performance, to simplify the analysis, it is 



assumed that the sensor’s transfer function has a unity 
value, implying a system’ s unity feedback. 

Since the digital computer can be programmed to 
perform any form of numerical computation, the digital 
compensator can realize any type of controller equation 
or algorithm known as the difference equation. Further, 
the digital control system shown in Figure 1 is modified, 
taking into account that the combination of the ADC, 
digital computer and DAC can be accurately represented 
by the combination of an ideal sampler, with a sampling 
period T , a digital filter D(z ), that solves the difference 
equation, and a zero-order hold (ZOH) [5, 6]. The 
complete modified model is shown in Figure 2. 



i i 




Figure 1. Basic Components of a Digital Control System 



I 1 




Figure 2. Modified Model of a Digital Control System 



For the purpose of the analysis and design, the modified 
system model of Figure 2 will be further simplified by 
showing only the sampler, the zero-order hold circuit and 
the plant, as shown in Figure 3. After the system’s 
analysis, a digital robust controller is designed in terms 
of its difference equations and is included in the structure 
of the system. 



For reaching the objectives of this research, a complex 
real-time digital tracking control system is considered. Its 
continuous part is a high order plant that is preliminary 
reduced and approximated to a third order system. The 
block diagram shown at Figure 3, is not including 
currently the digital filter that is assumed as D(z) = 1 
during the initial analysis of the system. 



R(s) + 



Sampler 



ZOH 



Plant 



"A E(s) / E(z) 


Gzoh(s) 


U(s) 


G P (s) 


C(s) 




















Figure 3. Further Simplified Digital Control System 



3. D-Partitioning Analysis in the Discrete- 
Time Domain 

The plant transfer function, as seen in Figure 3, is: 

G (s) = (1) 

p 0.0004 s 3 + 0.085i 2 + s 



The plant’ s gain K is considered as uncertain and variable 
system parameter. To determine the gain limitations and 
its marginal value, the method of the D-partitioning is 
applied in a similar manner like in its application for 
continuous systems [2, 6]. For the digital control system 
analysis, the D-partitioning is applied only to the 
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continuous plant, on condition that the sampling period T 
is taken into account. To achieve D-partitioning analysis 
accuracy, the sampling period T s should be at least 10 
times smaller than the plant’s minimum time constant 
T min , or T s < 0.\T min [1, 6]. For the discussed system, the 
sampling period is chosen as T s = 0.0001 sec, while the 
plant’s minimum time constant is T min = 0.005sec. 
Following the steps of the D-partitioning analysis, the 
characteristic equation of the continuous plant G P (s) is 
presented as: 

G(s) = 0.0004s 3 + 0.085s 2 +s + tf = 0 (2) 

From equation (2), the variable parameter K of the 
original system is determined as: 

K(s) = -0.0004s 3 - 0.085s 2 - s (3) 

The D-partitioning in terms of the parameter K , as seen 
in Figure 4, is achieved from equation (3) by following 
the procedure: 



» Hdlfb = feedback(Hdl,l) 

Transfer function: 

8. 82 8e -008 z A 2 + 3.513e-007 z + 8.735e-008 

z A 3 - 2.979 z A 2 + 2.958 z - 0.979 
Sampling time: 0.0001 
» step(Hdlfb) 



Marginal Step Response 




Figure 5. Marginal Step response for System Gain K = 213 



» K=tf([-0.0004 -0.085 -1 0],[0 1 ]) 
» nyquist(K) 



D-Partitioning in Terms of the Variable Gain 




Figure 4. D-Partitioning Analysis of the System in Terms of the 
Variable Gain K 



Taking into account that the frequency limits in the 
discrete-time domain are to = ± oo s /2 = ±2n/2T = ±31416 
rad/sec, the plotted D-partitioning curve in the ^f-plane is 
considered only within the range -31416 rad/sec < co < 
+31416 rad/sec [7, 8]. 

The D-partitioning, as seen in Figure 4, determines three 
regions on the ^f-plane: D(0), D(l) and D(2). Only D(0) 
is the region of stability, being always on the left-hand 
side of the curve [2, 9]. This implies that the system is 
stable only when its gain is within the range 0 < K < 213. 
At K = 213 the system becomes marginal. The result for 
the marginal system gain, obtained from the D- 
partitioning analysis can be confirmed in the discrete- 
time domain, Figure 5, in the following code: 

» Gpl=tf([0 213], [0.0004 0.085 1 0]) 

» Hdl = c2d(Gcl, 0.0001, 'zoh ') 



The optimal continuous plant gain ^ at a relative 
damping ratio 0.707 can be determined by taking into 
account the unity feedback of G P (s) and applying the 
following code [10]: 

» K=[0:0.001:10]; 

»for n=l:length(K) 

Gp_array(:,:,n)=tf([10000*K(n)],[4 850 10000 10000*K(n)]); 

End 

» I y,z]=damp(Gp_array ); 

» plot(K,z(l,:)) 



Relative Damping Ratio in Terms of the Gain K 




Figure 6. Determination of the Optimal Gain K 



The relative damping ratio ^ - 0.707 is achieved, at 
optimal continuous plant gain K. By interacting with the 
plot shown at Figure 6, the optimal gain is determined as 
K = 5.879. By applying the z-transform, the open-loop 
and closed-loop transfer functions of the system are 
determined in the discrete-time domain. The resulting 
time response to a step input signal is shown in Figure 7 : 



» Gp2=tf([0 5.879 ], [0.0004 0.085 1 0[) 
» Hd2 = c2d(Gp2, 0.0001, 'zoh ') 

» Hd2fb = feedback(Hd2,l) 



» step(Hd2fb )) 
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Figure 7. Step Response of the Optimal Digital System in Discrete-Time and its Zoomed-in Image 



The digital system performance is quite sensitive to the 
variation of the gain K. This can be demonstrated with 
system’ s time response for different gain values as shown 
in Figure 8: 



G p (s) 10000/j: 

l + G P (,s) "" 4 s 3 + 850s 2 + 10000.S + 10000 A" (4) 

58790 

4s 3 + 850s 2 + 10000s + 58790 




Figure 8. Step Responses of the Uncompensated Digital System in the 
Discrete-Time Domain 

4. Robust Controller Design 

The performance criterion “Integral of Time multiplied 
by the Absolute value of Error” (ITAE) is applied to 
achieve an optimal robust control. It is known that the 
ITAE reaches a minimum value for a relative damping 
ratio C, - 0.707 [3, 8, 9]. A solution, the system to meet 
the ITAE, is to implement a robust controller consisting 
of a series and a forward stage. The controller design 
employs a two-step compensation that operates in two 
degrees of freedom [10, 11, 12] and enforces a desired 
system damping, stability and time response. 



This condition corresponds to the desired closed-loop 
poles -6.05 ± j6.06. They are determined by the code as 
shown below: 

» Gcl=tf([ 58790], [4 850 10000 58790]) 

» damp(Gcl) 

Eigenvalue Damping Freq. (rad/s) 

-6.05e+000 + 6.06e+000i 7.07e-001 8.56e+000 

-6.05e+000 - 6.06e+000i 7.07e-001 8.56e+000 

-2.00e+002 1.00e+000 2.00e+002 

After approximating the system pole values to -6.1 ± j6.1, 
the two controller zeros can be placed at -6.1 ± j6.1. For 
physical realization of the controller’s series stage two 
remote controller poles, are added at s 1>2 = -1000, jO, so 
that their effect on the system performance is negligible 
[3, 14, 15]. Thus, the transfer function of the series robust 
controller stage is presented as: 

(5 + 6.1 + j6.1)0 + 6.1- y'6.1) 

G s ( 5 ) = = 

74.42(s + 1000) 2 
_ s 2 + 12.2s + 74.42 
74.42(s + 1000) 2 

However, to simplify the analysis, the remote poles s 1>2 = 
-1000, jO are not included into the transfer function of the 
open-loop plant: 



The design strategy for constructing the series stage G s (s) 
of the controller is such that its two zeros are placed near 
the desired closed-loop poles satisfying the optimal 
system condition. The optimal value of K corresponding 
to the relative damping ratio £ = 0,707 of the closed-loop 
system is already determined in the previous section as K 
= 5.879. Then, the unity feedback transfer function of the 
plant can be presented also as: 



g ol 0 s ) = G s ( s ) G p ( s ) = 



0.01344 (5 + 12.25 + 74.42) 
0.0004.S 3 + 0.085s 2 + s 



( 6 ) 



As a result, the unity feedback transfer function, 
including the series controller is: 
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0.01344 K(s 2 + 12.2 s + 74.42) 

G CL ( s ) — 

0.0004s 3 + 0.085s 2 + s + 0.01344 K(s 2 + 12.2s + 74.42) 



As seen from equation (7), the zeros of the closed-loop 
transfer function are in the vicinity of its poles. As a 
result, the closed-loop zeros will cancel the closed loop 
poles of the plant. To avoid this problem, a forward 
controller G F ( s) [10, 16] is added in series to the closed- 
loop transfer function G C l( s), as shown in Figure 9. 

The strategy for design of the forward stage G F ( s) is such 
that its poles should cancel the zeros of the closed-loop 
plant transfer function G C l( s). As a result, the total 
compensated part of the system Gj(s) will be dominated 
mainly by the two poles -6.1 ± j6.1, optimizing its 



( 7 ) 

performance. To improve on the step response rise time, a 
zero is added to the transfer function of the forward stage. 
Then, the transfer function of the forward controller stage 
is presented as: 

^ _ 74.42(0.15 + 1) 

Gp ( 5 ) — - 

The transfer function of the total compensated part, as 
seen from the block diagram in Figure 9, is presented as: 



^ ^ K(0 . 15 + 1) 

G t (s) — Gp(s^Gpp(s^ — (9) 

0.00045 3 + 0.085 5 2 + 5 + 0.01344^(5 2 + 12.25 + 74.42) 




Figure 9. Analogue Robust Controller Incorporated into the Control System 



The D-partitioning in terms of the variable parameter K 
is determined from the characteristic equation of Gj{s) 
and is plotted in Figure 10 following the course of action: 

- 0.0004 s 3 - 0.085s 2 - s 

K = 2 (10) 

0.0134s + 0.164s + 1 

»K = tf([-0.0004 -0.085 -1 0], [0.0134 0.164 1]) 

» nyquist(K) 

It is seen from Fig. 10 that the D-partitioning curve 
determines four regions of the ^-plane: D(0), D(l), D(2) 
and D(3). Only D(0) is the region of stability, being 
always on the left-hand side of the curve for a frequency 
variation from oo = -31416 to oo = +31416. It is seen that 
the system will be stable for any values of K > 0, since 
these values are within the region D(0). Also, the 
compensated system is examined for robustness in the 
Discrete-time domain by substituting different random 
values for the gain, K = 20, 50 and 200 in equation (9): 



D-Partitioning of the Robust Control System 




Figure 10. D-Partitioning of the Robust Control System in Terms of the 
Variable Gain K 



G Tl (s) 



20(0.15 + 1) 



20(0.15 + 1) 



Gj 9 (D 



0.0004 5 3 + 0.0855 2 + 5 + 0.27 (s 2 + 12.25 + 74.42) 0.00045 3 + 0.355 2 + 4.285 + 20 

50(0.15 + 5 ) 50(0.15 + 1) 

K = 50 0.0004 .j 3 + 0.085 .s 2 + s + 0.67 (s 1 + 12.2s + 74.42) 0.0004.? 3 + 0.76^ 2 + 9.2 s + 50 



Gy-O (S) 



200(0. li + 1) 



200(0. li + 1) 



( 11 ) 

( 12 ) 

(13) 



K ^ 200 0.0004s 3 + 0.085s 2 + s + 2.69(s 2 + 12.2s + 74.42) 0.0004s 3 + 2.77s 2 + 33.7s + 200 
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The system performance in the discrete-time domain, as 
shown in Figure 1 1, is obtained by the following code: 

GTl=tf([2 20], [0.0004 0.35 4.28 20]) 

GT2=tf([5 50], [0.0004 0.76 9.2 50]) 

GT3=tf([20 200], [0.0004 2.77 33.7200]) 

Hdl = c2d(GTl, 0.0001, ' zoh ') 

Hd2 = c2d(GT2, 0.0001, ' zoh ') 

Hd3 = c2d(GT3, 0.0001, ' zoh ') 

Hdlfb = feedback(Hdl,l) 

Hd2fb = feedback(Hd2,l) 

Hd3fb = feedback(Hd3,l ) 

step(Hdlfb,Hd2fb,Hd3fb) 

It is seen from Figure 1 1 that as a result of the applied 
two- stage robust controller, the digital control system 
becomes quite insensitive and robust to variation of the 
gain K. 

The step responses, shown in the discrete-time domain, 
for cases of K = 20, K = 50 and K = 200 differ 
insignificantly from each other. 



Step Responses of the Digital Robust Control System 




Figure 1 1 . Step Responses of the Digital Robust Control System in the 
Discrete-Time Domain 



damp(Hdfbl) 

Eigenvalue 


Magnitude 


Equiv. Damping 


Equiv. Ereq. (rad/s) 


9.99e-001 + 4.46e-004i 


9.99e-001 


8.10e-001 


7.61e+000 


9.99e-001 - 4.46e-004i 


9.99e-001 


8.10e-001 


7.61e+000 


9.17e-001 


9.17e-001 


1.00e+000 


8.63e+002 


damp(Hdfb2) 

Eigenvalue 


Magnitude 


Equiv. Damping 


Equiv. Ereq. (rad/s) 


9.99e-001 + 5.41e-004i 


9.99e-001 


7.46e-001 


8.14e+000 


9.99e-001 - 5.41e-004i 


9.99e-001 


7.46e-001 


8.14e+000 


8.28e-001 


8.28e-001 


1.00e+000 


1.89e+003 


damp(Hdfb3) 

Eigenvalue 


Magnitude 


Equiv. Damping 


Equiv. Freq. (rad/s) 


9.99e-001 + 5.93e-004i 


9.99e-001 


7.16e-001 


8.50e+000 


9.99e-001 - 5.93e-004i 


9.99e-001 


7.16e-001 


8.50e+000 


5.01e-001 


5.01e-001 


1.00e+000 


6.91e+003 



As seen from the results, the digital system has become 
quite robust. In spite of the considerable change of the 
gain, from K = 20, to K = 50 and further to K- 200, the 
system’s dominant conjugative poles and the system’s 
relative damping ratio change insignificantly. 

One of the approaches used to realize the series and 
forward robust controller stages in analogue systems is to 
implement RC networks and operational amplifiers [8, 17, 
18]. 

Alternatively, digital filters corresponding to the two 
robust stages can be designed and implemented in the 
system. 



Taking into account equation (5), the transfer function of 
the series digital robust stage can be determined by: 

» Gs=tf([l 12.2 74.42], [74.42 148840 74420000 ]) 

» Ds = c2d(Gs, 0.0001, 'zoh ') 

Transfer function: 

0.01344 z A 2 - 0.0268 z + 0.01336 

z A 2 - 1.81 z + 0.8187 
Sampling time: 0.0001 

To enable the implementation of a microcontroller in the 
mini-loop of the system, the transfer function of the series 
digital robust stage D s (z) is expressed in the discrete-time 
domain [19, 20, 21] by the following difference equation: 



y(lcT) = 0.01 34 4x(kT) - 0.0268 x[(k - l)r] + 0.0 1336 <x[(k - 2)T ] + 1 .8 ly [(k - l)r] - 0.8 1 87y[(/k - 2)T ] (14) 



Similarly, from equation (8), the transfer function of 
the forward digital filter stage is determined by: 

» Gf=tf([7.442 74.42], [1 12.2 74.42]) 

» Df= c2d(Gf 0.0001, 'zoh') 

Transfer function : 

0.0007441 z - 0.0007434 



To enable the implementation of a microcontroller in 
series with the mini-loop of the system the transfer 
function of the forward digital robust stage D F (z) is 
expressed in the discrete-time domain [19, 20, 21] by the 
following difference equation: 



z A 2- 1.999 z + 0.9988 
Sampling time: 0.0001 
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m(kT) = 0.0007441e[(& - l)T ] - 0.0007434e[(& - 2)r] + 1.999m[(& - \)T ] - 0.9988m[(& - 2 )t] (15) 



The block diagram of Figure 9 is modified to show the 
operation of the robust control system in the discrete- 
time domain, as seen in Figure 12. 

Each combination of a sampler, a digital filter D F (z) or 
D s (z ), solving a difference equation, plus a zero-order 
hold (ZOH), can be represented by a combination of 



the ADC, microcontroller and DAC, as shown in 
Figure 13. It is seen that the two microcontrollers, 
representing the series and the forward stages are 
incorporated into the control system. They are 
accordingly programmed to solve the difference 
equations (14) and (15) [6, 22]. 



R(s) + E(s) E(z) M(z) V(s) + X(s) X(z) Y(z) U(s) C(s) 




Figure 12. Two-Stage Digital Robust Controller Incorporated into the Control System 



ADC Microcontroller DAC ADC Microcontroller DAC 




Figure 13. Two-Stage Robust Microcontroller Incorporated into the Control System 



5. Conclusions 

Contributions of this research are the advancement of the 
D-partitioning method in its application for analysis of 
digital control systems and digital robust control design. 
The results of the research demonstrate that the D- 
partitioning analysis can be applied to the continuous 
plant in a digital control system, if the system’ s sampling 
period is at least 10 times smaller than the plant’s 
minimum time constant. If this condition is not met, the 
accuracy of the D-partitioning analysis is to some extent 
affected, but the method still can be applied with a 
reasonable effect. 

With the aid of the D-partitioning analysis, the marginal 
conditions of the digital system are determined and 
further the results are confirmed by the system’s step 
response in the discrete-time domain. 

The original digital control system is examined for 
robustness. It is established that the system is quite 
sensitive to any changes of the plant’s gain. Based on the 
ITAE criterion, a design strategy of a robust controller is 
suggested. 

An optimal robust controller is achieved by applying 
forward- series compensation with two degrees of 
freedom. The controller enforces a desired system 
damping and therefore a desired margin of stability. The 
system is analyzed with the robust controller and the 



results prove that it becomes quite insensitive to the gain 
variations. 

Following the robust controller design, two options of its 
implementation are suggested. One of the solutions is to 
construct the series and forward controller stages with the 
aid of RC circuits and operational amplifiers and applied 
them in the system as analogue continuous components. 
Alternatively, the transfer functions of the two controller 
stages are presented in the discrete-time domain by their 
difference equations. 

Two microcontrollers are incorporated into the control 
system as series and forward robust control stages and are 
programmed to solve their difference equations. Due to 
the advantages associate with the digital systems, the 
microcontroller robust design strategy is the better option 
and can be applied for any control system. 
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Abstract 

Wind turbines are complex, nonlinear systems sub- 
jected to large loads that reduce performance and 
lifetime. To ameliorate this, high performance control 
systems must be implemented. Plant characteristics 
such as flexible body dynamics and nonminimum 
phase reduce available control performance as their 
size increases. A comparative analysis of available 
control performance for a 600 kW and a 5 MW wind 
turbine is presented. Control limitations for each 
are identified, and available feedback is quantified. 
Dimensionless parameters are defined to provide 
response-based measures for comparison. Results 
of nonlinear simulations for both wind turbines 
are presented and compared. Aggressive control 
strategies are developed in an attempt to improve the 
performance of the larger wind turbine’s controllers. 

Keywords: Wind energy, Region 3 Control, 

Speed regulation, Vibration suppression, Single- 
input/ Two- output, Variable gain, Nyquist- stable con- 
trol 

1 Introduction 

Wind turbines are nonlinear systems, with character- 
istics that substantially restrict the available control 
system feedback: variable plant parameters, nonmini- 
mum phase zeros, adjacent pole pairs at low frequency, 
actuator delay and other limiting factors. Control per- 
formance, particularly in Region 3 (beyond rated wind 
speed), is critical, as the forces applied can severely 
reduce turbine lifetime. 

Wind turbine Region 3 rate regulation and load 
mitigation control has been extensively researched. 
Model-based gain selection strategies for PID Region 
3 compensators have been developed [1, 3]. Control 
of hybrid renewable energy systems are reported in 
[2]. Researchers at the National Renewable Energy 



Center (NREL) have designed and implemented ad- 
vanced control systems for rotor rate regulation and 
load mitigation applied to the Controls Advanced Re- 
search Turbine 2 (CART2), a two-bladed 600 kW 
Westinghouse wind turbine, for comparison to fixed 
gain low order controllers. Independent blade pitch 
control strategies to mitigate the effect of asymmetric 
wind variations across the rotor disk implemented on 
CART2 are reported in [4, 5]. Combined collective 
and individual blade pitch control for load mitigation 
is investigated in [6]. H in f control is used for speed 
and tower damping [7]. In [8], a Linear Quadratic 
Gaussian design with novel refinements is used to ac- 
count for CART2 actuator delay. State-space con- 
trol with disturbance accommodation is investigated 
in [9, 10, 11]. Full-state feedback is used for tower 
damping and torsion mode control for NREL’s three- 
bladed CART3 turbine in [12]. 

The US Department of Energy suggests wind 
power capacity will increase by roughly an order of 
magnitude (305 GW) by 2030 [13]. To meet this goal, 
deployed turbines must be substantially larger than 
what are deployed today. In [14], it is stated that 5 
MW turbines or larger are needed to meet the goals of 
the 2030 report. These turbines will have hub heights 
of 140 m and 153 m rotor diameters. Comparing 
this the 80 m hub heights and 77 m rotor diameters 
of currently deployed 1.5 MW turbines prompts two 
questions related to automatic control applications for 
these plants: what control performance is required for 
these 5 MW turbines, and what control performance 
is available? 

For a constant torque wind turbine, Region 3 con- 
trol has the primary goal of generator rate regula- 
tion. A low variance of this output is desired. Upper 
bounds on rate amplitude also apply subject to the 
limitations of the turbine power electronics (e.g. for 
CART2, a high speed shaft rate greater than 1860 
rpm is considered an overspeed event (rated speed is 
1800 rpm)). Ancillary Region 3 goals are related to 
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load mitigation. This includes using variable blade 
pitch to damp tower modes, and the use of generator 
torque to damp the shaft torsion mode. 

Control performance is directly proportional to 
feedback, which in turn is a function of the available 
control bandwidth. Some plant characteristics, like 
flexible body dynamics and nonminimum phase, have 
the effect of reducing available bandwidth. It is ex- 
pected that the flexible modes of proposed very large 
wind turbines (5 MW and larger) will be at lower 
frequencies than those of current turbines, and thus 
the available feedback is less. It is expected that the 
maximum blade angle rates will be lower than those 
of smaller turbines, also limiting the available control 
performance due to saturation. As a consequence, it 
is expected that the available control performance for 
the 5 MW turbines is less than that of current sys- 
tems. 

While effective, many of these control approaches 
make difficult the task of quantifying achievable per- 
formance. Modifying state or control penalties of an 
LQG controller effect changes in closed loop perfor- 
mance, but it is difficult to relate these changes to 
effects on performance or robustness in a quantifiable 
way. Adaptive control strategies can force system re- 
sponse to match that of a desired system, but there 
are potential robustness issues and it is not clear what 
the effect on disturbance rejection is. In this work, 
a comparative analysis is presented, quantifying the 
available control performance for the CART2 turbine 
and the 5 MW offshore wind turbine. Classical con- 
trol methods are employed for two purposes: these 
allow full use of plant frequency response data, and 
the control performance is easily quantifiable. The 
purpose of this paper is not to argue the relative mer- 
its of different control approaches as they relate to 
wind turbine applications, but to investigate trends in 
available control performance as these systems grow 
in size. 

The paper is organized as follows: the Region 
3 control architecture is presented, the available 
feedback for CART2 and the 5 MW wind turbine 
is quantified, the single input /two output control 
designs are detailed, and control simulation results 
using NREL’s Fatigue, Aerodynamics, Structure and 
Turbulence (FAST) models [15] for both systems 
are reported to illustrate a trend in control perfor- 
mance as wind turbine size increases. A discussion 
of aggressive control strategies for possible appli- 
cation to large wind turbines precedes the conclusions. 

Nomenclature. The rational function T(s) of the 
Laplace variable s is the loop transmission (alter- 
natively return ratio) of a feedback loop. \F(s)\ = 
|1 + T(s)| is the feedback. \F(s\ > 1, \F(s)\ < 1 and 
|T(s)| << 0 define negative , positive and negligible 
feedback, respectively. These definitions are rooted in 
the effect of feedback on the logarithmic response of 




Figure 1: Region 3 control block diagram. Torque 
Controller 1 is the Single Input/Single Output (SISO) 
system, R3 Pitch Controller is the Single Input/Two 
Output (SITO) system. 

the closed loop system to disturbances. For a low pass 
loop transmission, where |T(j^)| = 1 is the band- 
width (alternatively 0 dB crossover frequency ), and 
ujf where \T{joS)\ = Ao, Vcj < u;/, Ao a constant is 
the functional bandwidth. / = for uo in the set 
of nonnegative real numbers. Nonminimum phase is 
the phase lag not found using the Bode phase/gain 
relationship [16]. 

2 Control Architecture 

The goals of the control design are twofold. The pri- 
mary goal is turbine shaft rate regulation in Region 3 
(CART2 is a constant torque wind turbine, and for 
purposes of this comparison, the 5 MW wind tur- 
bine is assumed to be constant torque as well). Sec- 
ond, flexible body modes are targeted for feedback 
application for load mitigation. The control variables 
are collective blade pitch and generator torque. The 
outputs are high speed shaft rate and tower fore/aft 
acceleration. These are variables in two control sys- 
tems: a single-input/single-output controller (genera- 
tor torque input and high speed shaft rate output), 
and a single-input /two-output controller (collective 
blade pitch input, and high speed shaft rate and tower 
fore/aft acceleration are outputs). Figure 1 shows the 
block diagram for this controller. 

The generator and the high speed shaft are col- 
located, and the frequency responses for both the 
CART2 (figure 2) and the 5 MW turbine consist of 
alternating zero/pole pairs that facilitate high perfor- 
mance feedback design due to phase constraints. Low 
order feedback compensation for both turbines results 
in more than 20 dB attenuation of shaft vibrations at 
the modal frequency. This control has the added ben- 
efit of reducing the modulus of the high speed shaft 
response to collective blade pitch, making easier the 
task of gain stabilization of that mode in the shaft 
rate regulation design problem. 

The design of the collective blade pitch controller 
is more challenging, and is the focus of this compara- 
tive analysis. The plant is single-input /two-output 
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Generator/Rotor Rate 1 8 m/s 




Figure 2: CART2 Generator Torque to High Speed 
Shaft Rate Frequency Response at 18 m/s 



(SITO); the two outputs are high-speed shaft rate 
for the Region 3 regulation problem (primary), and 
tower fore/aft acceleration for the vibration suppres- 
sion problem (secondary). In contrast to the gener- 
ator/high speed shaft rate frequency response, the 
collective blade pitch responses present challenging 
design problems and it is this controller that distin- 
guishes the two turbines from the perspective of avail- 
able performance. 

The comparative analysis is presented thusly: 

1. Limitations to available feedback for the SITO 
collective blade pitch controller are identified 
and quantified for both the CART2 and 5 MW 
wind turbines. 

2. The same loop transmission shape is applied to 
both the CART2 and 5 MW high speed shaft 
rate control design that produces robust, high 
performance, linear control. 

3. Using 1 and 2, the available feedback over a cho- 
sen functional bandwidth for both turbine con- 
trol designs is calculated and compared. 

4. Controller designs for both systems are devel- 
oped, and the performance is compared using 
FAST simulations of both wind turbines. These 
results are compared to the trend indicated in 
3. 

It is assumed that adequate driveshaft vibration 
attenuation is achieved via simple generator torque 
control, and details of this design are not included in 
the comparative analysis. 



3 Comparison of Available Feed- 
back: CART2 vs 5 MW 

The performance difference between the collective 
blade pitch control performance of the 600 kW and 
5 MW wind turbines is quantified. The following con- 
ditions apply to both systems. 

1. The plant is two output: high speed shaft rate, 
and tower fore/aft acceleration. 

2. The control system is linear. 

3. The blade pitch rate is limited. 20 deg/s for 
CART2, 8 deg/s for the 5 MW turbine. 

4. The functional bandwidth, defined as the fre- 
quency to which the loop transmission has a 
constant magnitude Ao, is 0.1 rad/s. This allows 
more than one decade for the loop transmission 
roll off. 

5. The loop transmission roll off at frequencies 
higher than the functional bandwidth is -10 
dB/oct. This approach combines aggressive 
feedback delivery with good robustness (-150 de- 
gree phase after the break). 

6. If there are x octaves between the functional 
bandwidth and the 0 dB crossover frequency, 
then the loop transmission modulus over the 
functional bandwidth is 10(x+l) dB. This is ad- 
ditional 10 dB of feedback is extracted by careful 
loop shaping at the functional bandwidth break 
and near 0 dB crossover. 

The desired loop transmission function is given by 
the following function [16]. 

Tdes = e lo9{Ao)+ i lo9e{jf \ ( 1 ) 

where 9(jf) = . 1 , and / is the frequency nor- 

V 

malized with respect to functional bandwidth. Figure 
3 shows the modulus and argument of this function. 
Implemented controllers will be defined by rational 
functions with poles and zeros selected to approximate 
this shape. 

3.1 Available Feedback for the CART2 
Region 3 High Speed Shaft Rate 
Regulator 

Figures 4 and 5 show the CART2 high speed shaft rate 
and tower fore/aft acceleration responses to collective 
blade pitch at 18 m/s wind speed, found using the 
FAST model of this turbine. These plots are used to 
determine the feedback limitations for the high speed 
shaft rate controller, the primary controller for this 
application. 
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Figure 3: Desired loop transmission as a function of 
frequency normalized with respect to functional band- 
width. 



Collective Blade Pitch-Rotor Rate 18 m/s 




10 2 10 1 10° io 1 io 2 

Frequency (rad/sec) 

Figure 4: CART 2 Collective Blade Pitch to High 
Speed Shaft Rate Frequency Response at 18 m/s 



Collective Blade Pitch-Tower F/A Acceleration 18 m/s 




Figure 5: CART2 Collective Blade Pitch to Tower 
Fore/ Aft Acceleration Frequency Response at 18 m/s 



1. There is substantial nonminimum phase in the 
neighborhood of the first tower mode at 5.5 
rad/s of the shaft rate response (figure 4). This 
is the result of right-half plane zeros clustered 
near the first tower fore/aft and side-to-side 
modes at this frequency. The nonminimum 
phase delay at this frequency makes impossible 
the design of a shaft rate controller of higher 
bandwidth. 

2. There must be at least one octave of separation 
between the frequencies where the high speed 
shaft rate is negative and where the acceleration 
feedback is negative. 

3. The first blade edgewise mode at about 70 rad/s 
in the shaft rate response must be either gain 
or phase stabilized. The drivetrain mode at 
22 rad/s is reduced by the generator controller, 
however the collective blade pitch controller will 
stabilize this mode in case the generator loop is 
opened. 

Items 1 and 2 are the most substantial limitations. 
The nonminimum phase at 5.5 rad/s described in item 
1 requires gain stabilization of modes near this fre- 
quency. This limits the bandwidth of the high speed 
shaft rate controller to ~l-2 rad/s. Control of the 1st 
tower fore/aft mode (5.5 rad/s) requires acceleration 
feedback over the interval [4 6] rad/s, thus the restric- 
tion of item 2 limits the high speed rate feedback loop 
to 2 rad/s. The bandwidth limit caused by item 3 is 
higher than those of 1 and 2, and thus is not a critical 
limitation. 

The crossover frequency is ^1.5 rad/s. The inter- 
val from the functional bandwidth (0.1 rad/s) to the 
crossover frequency of 1.5 rad/s is 3.9 octaves wide. 
Assuming careful loop shaping at the break, the high 
speed shaft controller can deliver 10(3.9 + 1) = 49 dB 
of feedback over the functional bandwidth. 

3.2 Available Feedback for the 5 MW 
Region 3 High Speed Shaft Rate 
Regulator 

Figures 6 and 7 show the 5 MW high speed shaft rate 
and tower fore/aft acceleration responses to collective 
blade pitch at 18 m/s wind speed. These plots are 
used to list the feedback limitations for the high speed 
shaft rate controller. 

1. There is substantial nonminimum phase in the 
neighborhood of the first tower mode at 2 rad/s 
of the shaft rate response (figure 6). This is 
the result of right-half plane zeros clustered near 
the first tower fore/aft and side-to-side modes at 
this frequency. The nonminimum phase delay at 
this frequency makes impossible the design of a 
shaft rate controller of higher bandwidth. 
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Collective Blade Pitch-Rotor Rate 18 m/s 




trol of the 1st tower fore/aft mode (2 rad/s) requires 
acceleration feedback over the interval [1 4] rad/s, 
thus the restriction of item 2 limits the high speed rate 
feedback loop to 0.5 rad/s. The driveshaft and edge- 
wise modes are sufficiently higher frequency than the 
limits of items 1 and 2 that item 3 is not a limitation. 

The crossover frequency is ^0.5 rad/s. The in- 
terval from the functional bandwidth (0.1 rad/s) to 
the crossover frequency of 2 rad/s is 2.3 octaves wide. 
Assuming careful loop shaping at the break, the high 
speed shaft controller can deliver 10(2.3 + 1) = 33 dB 
of feedback over the functional bandwidth. 

The more substantial limitations of the 5 MW tur- 
bine result in a factor of approximately six times less 
rate regulation performance than the CART2, theo- 
retically. 



Figure 6: 5 MW Collective Blade Pitch to High Speed 
Shaft Rate Frequency Response at 18 m/s 



Collective Blade Pitch-Rotor Rate 18 m/s 




Figure 7: 5 MW Collective Blade Pitch to Tower 
Fore/ Aft Acceleration Frequency Response at 18 m/s 



4 SITO Designs for CART2 and 
5 MW 

Single-input /single-output (SISO) generator con- 
trollers are designed for each turbine to damp the 
shaft torsion mode. These are simple controllers and 
are not described in this work. The collective blade 
pitch SITO designs are subject to the bandwidth lim- 
itations quantified in the previous sections, and clas- 
sical loop shaping methods are employed to extract 
maximum feedback over the functional bandwidth. It 
is noted that although the controllers are designed us- 
ing linearized dynamic models at a fixed wind speed, 
they will be tested in turbulent wind conditions over a 
broad interval of Region 3 speeds. The compensators 
for the CART2 SITO collective blade pitch controller 
(high speed shaft rate G rate (C 2 ) i s ) and structure con- 
trol G struct{C 2 ){s)) are 

r / \ 0.1(s 2 + a + 0-25) . 

rate(C2) V S 4 + 2 0.1s 3 + 102s 2 + 10.2s + 1 U 



2. There must be at least one octave of separation 
between the frequencies where the high speed 
shaft rate is negative. 

3. The first blade edgewise mode at about 25 rad/s 

in the shaft rate response must be stabilized (ei- 

ther gain or phase). The drivetrain mode at 10 
rad/s is attenuated by the generator controller, 
however the collective blade pitch controller will 
stabilize this mode in case the generator loop is 
opened. 



G s truct(C2) (^) \ a' (^) 

v ' S + D 

The compensators for the 5 MW SITO collective 
blade pitch controller (high speed shaft rate G r {s) and 
structure control G struct ^(s)) are 

0.09(s 3 + 3.5s 2 + 3.5s + 1) 

_ s 5 + 60.1s 4 + 531s 3 + 1303s 2 + 130.3s + 12.5 

(4) 



As is the case for CART2, items 1 and 2 are the 
most substantial limitations. While this response does 
not have the sharp poles that CART2 has around 2 
rad/s, the nonminimum phase at this frequency de- 
scribed in item 1 requires gain stabilization of modes 
near this frequency. This limits the bandwidth of the 
high speed shaft rate controller to ^0.5 - 1 rad/s. Con- 



0.9525(s 2 + 19s + 84) 
G st ruct(5){s) ~ s 3 + 58 s 2 +416s + 800 - 

The loop transmissions for these controllers are 
shown in figure 8. The first 0 dB crossover frequen- 
cies and functional bandwidth feedback are consistent 
with the projected maximum values determined in 
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Type 0 Loop Transmission 




Frequency (rad/sec) 



SITO Loop Transmission 




Figure 8 : Region 3 Collective Blade Pitch SITO Loop 
Transmissions for the CART 2 (blue) and 5 MW (red) 
wind turbines at 18 m/s 



Figure 9: Type 1 Region 3 Collective Blade Pitch 
SITO Loop Transmissions for the CART 2 (red) and 5 
MW (blue) wind turbines at 18 m/s 



section 3. The relative stability of both controllers 
are comparable; there is approximately 6 dB of gain 
stability margin and 30 degrees of phase stability mar- 
gin for both controllers. 

5 Control Performance with 
Turbulent Wind 

The analysis of the previous section is tested in sim- 
ulation using FAST models of these two turbines in 
turbulent wind environments, where wind speeds vary 
between 10 and 27 m/s. The loop transmission func- 
tions are slightly modified to include origin poles for 
asymptotic rejection of step wind disturbances (see 
figure 9). It is noted that while there no longer is a 
functional bandwidth for these controllers, the differ- 
ence in loop transmission modulus at 0.1 rad/s is the 
same as what is reported in section 3. 

To quantify the relative performance of the SITO 
controllers, two dimensionless parameters are defined. 
The high speed shaft rate standard deviation to rated 
speed ratio (SDRS) is 

SDRS = — — — , ( 6 ) 

^ hss(rated ) 

where ahss and u hs s(rated) are the high speed shaft 
rate standard deviation and the Region 3 rated high 
speed shaft rate, respectively, in the same units. As 
the turbines in this study are constant torque, SDRS 
provides a measure of power variation in Region 3. 

The secondary goal of the collective blade pitch 
SITO controller is tower fore/ aft acceleration attenu- 
ation. The tower fore/aft acceleration ratio (TFAAR) 
is 



Table 1 : SITO Control Performance in Turbulent 
Wind 





CART 2 


5 MW 


Rated Shaft Speed (rpm) 


1800 


1174 


HSS rate std (rpm) 


31.9 


83.46 


Tower Acc std (ol) (m/s/s) 


0.3549 


0.2988 


Tower Acc std (cl) (m/s/s) 


0.2193 


0.1719 



TFAAR =^l t ( 7 ) 

^a(ol) 

where <r a ( c q is the tower fore/aft acceleration standard 
deviation with the SITO controller closed, and <r a ( 0 n 
is the tower fore/aft acceleration standard deviation 
with high speed shaft rate feedback only (referred to 
as the baseline). 

FAST simulations are performed for both turbines 
over 10 minutes in turbulent wind. Table 1 shows the 
rated high speed shaft rate, the high speed shaft rate 
standard deviation, and tower fore/aft acceleration in 
open and closed loop conditions for both turbines with 
the SITO controllers with loop transmissions shown in 
figure 9. 

The TFAAR for the the CART 2 and 5 MW tur- 
bines are 0.62 and 0.58, respectively. The performance 
is comparable, which is consistent with the fact that 
both controllers apply close to the same feedback at 
the tower first mode frequency. The SDRS for the 
the CART 2 and 5 MW turbines are 0.018 and 0.071, 
respectively. The 5 MW SDRS is four times greater 
than that of CART 2 . This is consequence of the 5 
MW wind turbine plant dynamics restricting available 
feedback to greater effect than for the CART 2 turbine. 
As such, alternative methods, perhaps including pas- 
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sive damping, may be indicated for structural quieting 
of the 5 MW turbine. 

6 Aggressive Feedback Strate- 
gies for 5 MW SITO Designs 

The 5 MW wind turbine plant dynamics severely re- 
strict available feedback, which is quantified in section 
3. These results are consistent with other simulations 
with step wind changes, and more severe turbulent 
wind [17]. Two advanced control strategies are inves- 
tigated in this section in an attempt to improve the 
performance despite these limitations. 

6.1 HSS Rate Error Variable Gain 
SITO Controller 

The requisite frequency separation for the SITO col- 
lective blade pitch controller (item 2 from 3.2) is a sub- 
stantial limitation. As rotor rate regulation is the pri- 
mary function of the collective blade pitch controller 
in Region 3, a variable gain strategy is proposed to 
apply more rate feedback when error is large. This 
comes at the expense of lowering tower acceleration 
feedback to allow increased rate feedback bandwidth. 

The variable gain algorithm is a function of high 
speed shaft rate (c Uhss) = 1174 rpm and gear ra- 
tio r g = 97. The output is a gain applied to 
the nominal tower acceleration feedback compensator, 
kstruct(hv)(whss) and the nominal high speed shaft 
rate feedback compensator, k rate ( 5v \(ujhss)- Variable 
gain k struct ($ v }(ujhss) is defined in table 2. The con- 
stants l O h ss (i),i = 1,2, 3, 4 are threshold values. If 
the high speed shaft rate is in the interval bounded 
by ^hss( 2 ) and ^hss( 3 ) 5 the gain ^strnct( 5 v) is 

unity, and the collective blade pitch controller is the 
SITO controller previously described. If the shaft 
speed is outside this interval, k struct ^(ujhss) is re- 
duced linearly until the shaft speed is outside the in- 
terval bounded by 1 ) and ^hss( 4 ) where this gain 
is fixed at 0.05. Given this gain, the rate regulator 
gain is ^yate( 5 i>) (^/iss) ^k S f ruc i(j )V ^{uj] ljSS ^ T 5. 

For the 5 MW wind turbine, the threshold val- 
ues are c o hss (i) are 1039, 1125, 1222, 1309 rpm for i = 
1,2,3, 4. Figure 10 shows the loop transmission of the 
variable gain controller with different high speed shaft 
rate errors. The loop transmission associated with the 
largest rate error (light blue) has negligible feedback 
at the first tower mode frequency (approximately 2 
rad/s), and has about one octave more bandwidth (1 
rad/s) than the nominal SITO controller (0.4 rad/s). 
It is noted that the gain changes must be slower than 
the dynamics of the plant. A low pass filter at the 
output of the the variable gain system can be used to 
adequately slow the gain change [18]. 

Figures 11 and 12 show the high speed shaft rate 
for the 5 MW turbine with high speed shaft rate feed- 



Table 2: HSS Rate Error Variable Gain Algorithm 



^ hss 


k s truct(5v) {^hss^j 


M hss £ ^hss( 1) 


0.05 


^hss( 1) ^ ^JJiss 2l ^hss( 2) 


~^~{^hss 


^hss( 2) ^ ^ hss ^ ^hss( 3) 


1.0 


^hss( 3) ^ ^ hss ^ ^hss( 4) 


— ^{Uhss — U hss ( 3)) + 1 


^ hss ^ ^hss( 4) 


0.05 



Bode Diagram 




Frequency (rad/sec) 

Figure 10: Variable Gain 5 MW SITO Loop Trans- 
mission. Dark blue is the standard SITO loop trans- 
mission. Green, red and light blue loop transmission 
are in conditions of increasing high speed shaft rate 
error. 



5 MW HSS Rate 




Figure 11: 5 MW high speed shaft rate for base- 
line (blue), fixed gain SITO (green) and variable gain 
SITO (red) collective blade pitch control. 
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J i i i i i i i W i i i 

550 555 560 565 570 575 580 585 590 595 600 

time (s) 

Figure 12: 5MW high speed shaft rate for baseline, 
fixed gain and variable gain collective blade pitch con- 
trol. 1250 rpm is the variable gain threshold value. 

back only (baseline), fixed gain SITO and variable 
gain SITO control in turbulent wind conditions. It 
is evident that variable gain system has superior per- 
formance compared to the fixed gain SITO controller. 
The SDRS for the variable gain controller is 0.0575, 
a 19% improvement over the fixed SITO controller. 
Figure 13 shows the tower fore/ aft acceleration for 
the same controllers. As the acceleration feedback is 
reduced as a function of rate error, the variable gain 
controller has inferior performance to the fixed gain 
SITO system for tower load mitigation, but clearly 
superior to the baseline controller. The TFAAR for 
the variable gain controller is 0.67, a 16% reduction 
in performance compared to the fixed gain controller. 
This relative performance is a consequence of the vari- 
able gain strategy reducing the tower mode feedback 
when shaft rate error exceeds the threshold. 

6.2 Nyquist-stable Controller 

Another alternative to improve shaft rate regulation 
subject to the bandwidth limitations quantified in sec- 
tion 3 is to have a steeper roll off to the crossover fre- 
quency (i.e. roll off steeper than -10 dB/oct). This 
provides greater negative feedback over the functional 
bandwidth. A Nyquist-stable control system is one 
that is stable in closed loop, and has a stable loop 
transmission function that crosses the negative real 
axis of the Nyquist plot outside the unit circle. This 
describes a loop transmission function that over an in- 
terval of frequencies has a modulus slope steeper than 
-12 dB/oct. Nyquist-stable controllers have been im- 
plemented on aerospace plants with excellent results 
[19, 20]. The tradeoff for this increased feedback is 
a susceptibility to oscillation in actuator saturation. 
Nonlinear dynamic compensation can be employed 
with a Nyquist-stable system to satisfy the conditions 



i i i i i i 

550 555 560 565 570 575 

time (s) 

Figure 13: 5MW tower fore/aft acceleration for base- 
line, fixed gain and variable gain collective blade pitch 
control. 

of absolute stability [20, 21]. 

To test the efficacy of a Nyquist-stable design 
for the 5 MW wind turbine, a SISO collective blade 
pitch/high speed shaft rate controller with this prop- 
erty is designed and tested in a turbulent wind envi- 
ronment. This controller has the advantage of the first 
tower mode control being abandoned (item 2 from sec- 
tion 3.2 is no longer a factor), so the bandwidth can be 
increased to 1 rad/s. Figure 14 show the Nichols plots 
without M-rings for the high speed shaft rate feedback 
only loop transmission (compensator is G rate (C 2 ) de- 
fined in section 2 4, the structure loop is opened) and 
the Nyquist-stable controller. 

The Nyquist-stable controller is extremely aggres- 
sive. In addition to the steep roll off, the stability 
margins are smaller than those of the SITO rate reg- 
ulation controller, and the torsion mode is phase sta- 
bilized as opposed to gain stabilized. These features, 
in addition to the extra octave of bandwidth, give the 
Nyquist-stable controller about 40 dB more feedback 
at 0.1 rad/s than the SITO rate regulator. 

The aggressive features of the Nyquist-stable con- 
trol system result in poor robustness. This is a par- 
ticulary bad characteristic for wind turbine applica- 
tions, as the plant dynamics vary with wind speed. 
Figure 15 shows the Nyquist-stable loop transmission 
for four different wind speeds in Region 3. While re- 
taining stability at wind speeds higher than the design 
speed (18 m/s), there are critical point encirclements 
at both the initial crossover and the phase stabilized 
torsion mode. FAST simulations with this controller 
with turbulent wind frequently show violent oscilla- 
tions within one minute of operation. 

The Nyquist-stable robustness problem can be ad- 
dressed via the application of gain-scheduling, how- 
ever this is an arduous task. In addition to schedul- 
ing the linear controller, requisite nonlinear dynamic 
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Loop Transmission 20 m/s 




Figure 14: 5MW Nyquist-stable Loop Transmission 
Nichols Chart. 



Nichols Chart 




Figure 15: 5MW Loop Transmission Nichols Chart for 
Multiple Wind Speeds. 



compensation for stability retention in actuator sat- 
uration would also have to be gain scheduled. The 
result would be a very complicated, perhaps fragile 
control system. In addition, the increase in perfor- 
mance at low frequency is offset by increased positive 
feedback near crossover, which might result in inferior 
SDRS performance. The advantages of Nyquist-stable 
control systems are realized when there is substantial 
feedback bandwidth, a characteristic that wind tur- 
bines do not have (item 1 from section 3.2). Thus, 
this type of aggressive, nonlinear control (Nyquist- 
stable control with nonlinear dynamic compensation) 
is ill-suited for wind turbines. 



7 Results 

A comparative analysis of available feedback for 600 
kW and 5 MW wind turbines for multiple output, Re- 
gion 3 control has been presented. Subject to feed- 
back limitations that are driven by plant dynamics, 
it is determined that for the same loop transmission 
shape, the 5 MW high speed shaft rate regulator 
controller delivers 16 dB less feedback over a func- 
tional bandwidth of 0.1 rad/s. To present an alter- 
native method of comparison, two performance mea- 
sures are introduced using stochastic parameters and 
system responses: the standard deviation rated speed 
(SDRS) ratio indicates high speed shaft rate stan- 
dard deviation normalized to rated shaft rate, and 
the tower fore/aft acceleration ratio (TFAAR), the 
ratio of closed loop to open loop tower fore/aft accel- 
eration standard deviation. FAST simulation results 
show that while the CART2 and 5 MW controllers de- 
liver comparable tower acceleration attenuation, the 
SDRS of the 5 MW controller is 4 times greater than 
that of CART2. 

In an attempt to close the gap in performance, ag- 
gressive control strategies are applied to the 5 MW 
control problem. A rate error variable gain algorithm 
is applied to the 5 MW SITO controller, effectively 
trading away acceleration feedback for increased rate 
feedback to improve the SDRS. Using this approach, 
a 19% decrease in SDRS is achieved at a cost of a 
16% increase in TFAAR. A Nyquist-stable controller 
is designed without acceleration feedback represent- 
ing a more aggressive approach. While the feedback 
is substantially increased subject to the limitations 
identified, plant parameter variations due to changing 
wind speed drive this controller unstable. Extremely 
carefully designed gain scheduling and nonlinear dy- 
namic compensation are required, and it appears that 
this system with its intrinsic robustness limitations 
would not be feasible for this application. 

8 Conclusion 

Frequency separated, two output control provides 
comparatively less rotor rate regulation for a 5 MW 
turbine compared to a 600 kW turbine. It is unknown 
what the consequence of this comparative lack of con- 
trol performance using this design approach for the 5 
MW wind turbines is. For instance, a constant torque 
turbine (e.g. CART2) requires high performance ro- 
tor rate regulation in Region 3. A rotor rate of 3% 
above rated speed is considered an overspeed event for 
this turbine. Aggressive feedback control is required 
to prevent overspeed events for this turbine in very 
windy conditions. As less feedback is available for the 
5 MW turbine, it is unlikely that it can operate in 
very harsh conditions with comparable performance, 
as this analysis suggests. However, if the 5 MW tur- 
bines are operated in a constant power mode, then 
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the torque control can reduce generator torque as the 
turbine speed increases, reducing the stringent rate 
control requirement. The 5 MW turbines will be de- 
ployed offshore, so the increased rotor noise at these 
high rates is not an issue. This perhaps would allow 
more feedback to be applied for structural load miti- 
gation. 
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